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PREFACE 


Although  the  MAGIC  code^  is  designed  to  be  user  friendly,  our 

experience  has  been  that  a  new  user  will  typically  require  several  years 
to  become  proficient  in  its  use.  This  simply  reflects  that  performing 

particle-in-cell  (PIC)  code  simulations  on  complex  physical  problems  is  an 
art.  Specifically,  it  requires  an  intimate  knowledge  of  mathematical 
models  and  their  application  in  particular  physical  regimes,  it  requires 
judgment  in  the  trade-offs  which  are  necessary  to  perform  practical 

simulations,  and  it  requires  the  ability  to  interpret  results  and  to 
recognize  failures  due  to  modeling  errors,  algorithm  incompatibilities, 

numerical  instabilities,  et  cetera,  as  well  as  outright  code  errors. 

The  purpose  of  the  canonical  problems  presented  in  this  report 

is  to  accelerate  the  learning  process.  Specifically,  we  have  developed 

twenty-two  simulation  problems  which  relate  to  pulsed  power  transport. 
These  problems  are  divided  into  the  following  three  sets: 

(1)  Electromagnetic  Propagation  -  simulating  a  coaxial  line  in 
cold  test. 

(2)  Particle  Kinematics  -  simulating  gyrokinetics,  space-charge 
limiting,  and  magnetic  insulation. 

(3)  Power  Transport  -  simulating  self-consistent  pulsed  power 
transport  processes. 

These  three  sets  of  problems  are  described  in  this  report  in  Sections  1, 

2,  and  3,  respectively.  The  problems  are  intended  to  be  solved 

sequentially,  since  many  of  the  effects  are  interrelated  and  the  problems 
gradually  increase  in  level  of  difficulty. 


t  R.  Goplen.  R.  E.  Clark,  3.  McDonald,  and  W.  M.  Bollen.  "User's  Manual 
for  MAGIC/  Version-September  1993,"  Mission  Research  Corporation 
Report,  MRC/WDC-R-069,  September  1993. 


The  individual  problems  were  drawn  from  several  sources; 
however,  most  were  developed  specifically  for  the  present  application. 
Similarly,  the  authors  are  responsible  for  most  of  the  analytical 
solutions*  exceptions,  such  as  the  8  x  ®  drift  and  field  emission,  are 
denoted  by  references. 


For  each  of  the  problems,  we  have  included  selected  simulation 
results,  and  have  also  included  complete  copies  of  the  input  data  used. 
All  simulations  described  in  this  report  were  performed  with  one  version 


( 3une  1983)  of  MAGIC. 


For  future  versions  of  the  code,  certain  data 


commands  may  chanqe  in  detail;  however,  it  should  be  possible  at  any  time 
to  obtain  equivalent  results. 


Clearly  it  is  not  possible  to  treat  computational  pulsed  power 
transport  in  a  comprehensive  manner  by  performinq  a  few  simulation 
problems.  However,  it  should  be  possible  for  a  new  user  to  work  all  of 
these  problems  in  a  relatively  short  time  period  (say,  three  months);  this 
was  considered  to  be  an  important  constraint.  A  second  constraint  was 
imposed  by  available  computer  time:  the  exercises  would  need  to  be 
affordable.  Thus,  a  maximum  limit  of  five  minutes  CPU  time  on  the  CRAY  1 


was  set  for  any  simulation.  This  effectively  precluded  realistic 
self-consistent  problems.  Because  of  such  constraints,  this  experience 
can,  at  best,  offer  some  exposure  to  the  major  ideas. 


Finally,  the  authors  wish  to  thank  Games  Benford  of  Physics 
International  Company  for  suqgestlng  the  canonical  problem  approach,  and 
John  Brandenburg  of  NFC  for  a-sistance  with  the  analytical  solutions  in 


Problems  2.6  and  3.6. 
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SECTION  1 

ELECTROMAGNETIC  FIELDS 


The  first  problem  set  consists  of  purely  electromagnetic  field 
problems.  Since  Maxwell's  equations  in  the  absence  of  space  charge  are 
linear,  this  greatly  simplifies  the  interpretation  of  results  and  provides 
an  ideal  means  of  introducing  PIC  simulation  concepts.  Similarly,  Section 
2  will  begin  by  studying  the  motion  of  single  particles  in  specified  elec¬ 
tromagnetic  fields,  which  is  also  easy  to  check  analytically.  The  more 
complex  nonlinear  processes  involving  space  charge  will  be  deferred  until 
Section  3. 

The  electromagnetic  field  problems  are  all  based  upon  a  coaxial 
transmission  line  geometry.  They  are  designed  to  offer  some  exposure  to 
the  variety  of  algorithms,  models,  and  boundary  conditions  available  in 
the  code. 


•A? 

A* 

•  V  v 


COAXIAL  LINE. 


1.1.1  Problem  Description. 


Consider  the  section  of  coaxial  transmission  line  illustrated  in 
Figure  1.1.  The  inner  and  outer  radii  and  length  of  the  section  are 


defined  to  be 


r ^  =  0.10  m 


rn  =  0.20  m 


(1.1) 


L  =  0 . 60  m  . 

It  is  assumed  that  the  walls  are  perfectly  conducting  and  that  no 

free-charge  (e.g.,  from  field  emission)  can  be  created. 

For  the  simulation,  we  Introduce  at  the  inlet  an  incident 
voltage  pulse  representing  a  transverse  electromagnetic  (TEM)  wave 
traveling  from  left  to  right.  The  pulse  is  defined  to  reach  a  peak 

voltage  $+p  of  10 6  V  during  a  linear  rise-time  x  of  1  nsec  and  to  be 
constant  thereafter.  Thus,  the  incident  wave  equation  is 


0+(t)  = 


4>+(t/x)  ,  0  <  t  <  x 

P 


,  x  <  t  < 


(1.2) 


In  general,  the  wave  at  a  boundary  consists  of  both  incident  and  scattered 


parts,  or 


<t>  =  <t>+  + 


(1.3) 


However,  for  this  simulation,  we  want  the  scattered  wave  to  vanish 
( i*  =  0) ,  so  that  the  inlet  wave  will  consist  purely  of  incident  wave.  At 
the  outlet,  the  exact  opposite  is  to  be  the  case.  We  want  the  incident 


wave  to  vanish  (  i>+  =  0) ,  so  the  outlet  wave  will  consist  purely  of  scat¬ 
tered  wave.  This  Is  the  simplest  choice  of  boundary  conditions  which  will 
allow  us  to  watch  an  electromagnetic  wave  propagating  through  the  coaxial 
line. 


During  the  simulation,  we  wish  to  measure  the  voltage  <t>(t)  and 
the  current  I(t)  as  a  function  of  time  at  three  axial  locations:  z  =  0, 
L/2,  L.  We  also  desire  to  examine  the  radial  dependence  of  the  fields 
Er,  Ez,  and  Bg  at  the  same  three  axial  locations  at  the  end  of  the 
simulation,  which  is  to  cover  a  time  span  of  ten  nsec. 

1.1.2  Suggested  Approach. 

All  data  sequences  and  conventions  are  described  in  the  Septem¬ 
ber  1983  version  of  the  MAGIC  User's  Manual^.  Input  data  is  entered  in 
the  form  of  commands,  each  of  which  contains  a  key  word  followed  by  argu¬ 
ments  and  parameters  and  is  terminated  by  a  slash.  The  input  data  for 
this  particular  problem  are  given  in  Table  1.1.  For  this  first  problem, 
we  shall  briefly  discuss  all  of  the  input  data  in  the  order  of  appearance 
in  Table  1.1.  (Note  that  the  order  of  commands  Is  generally  arbitrary.) 
In  subsequent  problems,  only  important  variations  or  newly  introduced 
commands  will  be  discussed.  The  reader  must  refer  to  the  Manual  for  a 
detailed  description  of  command  parameters. 

The  TITLE  command  (see  Table  1.1)  provides  a  unique  problem 
Identification  on  all  output,  while  the  COMMENT  command  allows  the  user  to 
record  arbitrary  comments.  The  cylindrical  (z,r)  coordinate  system  is 
specified  by  the  SYSTEM  command,  which  Identifies  the  spatial  coordinates 
(xi  =  z,  X2  =  r).  To  resolve  the  spatial  coordinates,  we  shall  use  a  uni¬ 
form  grid  with  6r  =  6Z  =  0.01  m.  The  spatial  grid  must  be  sufficient¬ 
ly  fine  to  resolve  the  wave  front  and  associated  wavelengths.  The  z-axis 

t  B.  Goplen,  R.  E.  Clark,  3.  McDonald,  and  W.  M.  Bollen,  "User's  Manual 
for  MAGIC/Version-September  1983,"  Mission  Research  Corporation  Report, 
MRC/WDC-R-068 ,  September  1983. 
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is  specified  by  a  X1GRID  command  and  the  r-axis  by  a  X2GRID  command.  Note 
that  the  first  cell  (outer  boundary)  in  MAGIC  is  unusable  -  thus  the 
simulation  boundary  must  begin  at  spatial  index  2  rather  than  1.  We  use 
the  FIELDS  command  to  specify  a  centered-dif ference  field  algorithm  with 
500  time  steps  of  2  x  10“  sec  each,  for  a  total  duration  of  10'  sec. 
The  COURANT  command  specifies  a  test  for  compliance  with  the  Courant 
stability  criterion. 

The  CONDUCTOR  command  is  used  to  specify  the  conducting 
surfaces.  The  two  conducting  segments  are  arbitrarily  named  CATHODE  and 
ANODE.  Note  the  definition  of  the  surface  normal  flags  -  the  convention 
is  that  the  unit  vector  must  lie  normal  to  the  surface  and  point  to  the 
interior  of  the  simulation. 

The  VOLTAGE  command  is  used  to  specify  the  incident  voltage 
pulse  at  the  inlet.  In  particular,  the  option  which  zeros  out  scattered 
waves  at  the  inlet  (4)'  =  0)  is  specified  by  setting  IV0L(1)  to  zero.  (For 
this  problem,  there  should  be  no  physically  meaningful  scattered  wave.) 
Both  the  temporal  and  radial  distributions  of  the  incident  pulse  must  be 
specified,  and  these  are  simply  given  the  arbitrary  names,  TEMPORAL  and 
RADIAL,  respectively.  All  functions  in  MAGIC  are  entered  as  input  data 
using  FUNCTION  commands.  Thus,  the  temporal  function  in  Equation  (1.2)  is 
entered  in  the  form  of  numerical  data. 

The  radial  distribution  is  determined  as  follows.  It  can  be 
shown  that  this  system  of  two  separate  conductors  will  support  a  TEM 
waveA  This  wave  propagates  axially  at  the  speed  of  light,  has  only  two 
transverse  components  (Er  and  Bq),  and  has  a  transverse  field 
distribution  given  by  the  static  solution, 


y.'i 


8* 


1 


A 


»/  A 


AA 


t  3.  D.  Jackson,  Classical  E lectrodvnamics ,  3ohn  Wiley  and  Sons,  New  York, 
1967  (p.  243). 
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(1.4) 


.V.V.TT.V.  . 


V**  =  0  , 


which  in  cylindrical  coordinates  yields  the  results, 


E  =  -  - - - 

r  r An(r  /r . ) 


(1.5) 


In  this  equation,  <t>  is  the  anode  voltage.  For  the  radial  distribution, 
MAGIC  requires  only  the  relative  dependence  (r_1),  since  the  voltage 
normalization  is  specified  by  the  temporal  function.  The  radial  function 
is  entered  using  a  power  term  option. 

Finally,  the  LOOKBACK  command  is  used  to  specify  the  outlet 
boundary  condition  and  will  allow  outgoing  waves  to  escape  at  the  outlet. 
The  phase  velocity  factor  is  assumed  to  be  unity. 

To  measure  simulation  results^",  we  specify  time  histories  of 
the  radial  electric  field  integral  (voltage)  at  inlet,  midpoint,  and  out¬ 
let,  using  the  OBSERVE  command.  Six  magnetic  field  measurements  (cathode 
and  anode  at  the  same  axial  locations)  are  also  specified.  To  record  the 
spatial  variation  of  these  same  fields,  range  plots  are  specified  using 
the  RANGE  command.  The  OUTPUT  command  is  used  to  specify  output  on  system 
graphics  rather  than  the  Internal  MAGIC  line  printer  graphics. 

The  START  command  is  used  to  initiate  the  simulation,  and  the 
STOP  command  will  terminate  execution  after  completing  the  simulation. 

1.1.3  Analytical  Solution. 

Axial  propagation  in  a  coaxial  cable  can  be  considered  in  terms 
of  an  equivalent  transmission  Line.  Transmission  line  equations  for  TEM 


t  The  input  data  provided  for  these  problems  will  generally  produce 
copious  output;  for  reasons  of  space,  only  selected  simulation  results 
are  included  in  this  report. 
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propagation  involve  axial  and  temporal  derivatives  of  anode  voltage,  <fc, 
and  anode  current,  I.  The  equations  are 


a  <t>  =  l'5  i 
z  t 


a  i  =  C3  i 
z  t 


(1.6) 


For  a  coaxial  line,  the  capacitance  and  inductance  per  unit  length  are 


in(r  /r  ) 
o  I 


u  in(r  /r, ) 


(1.7) 


where  e0  and  Hq  are  the  free-space  permitivity  and  permeability, 
respectively.  Division  of  the  two  relations  in  Equation  (1.6)  yields 


4-  =  IZ  , 


where  the  line  impedance,  Z,  is 


Z  =  (L'/C')1/2  =  Z  Mr  /r.) 

o  o  i 


Z  =  _1  (jx  /e  )  1/2  =  59.96  ohms 
0  2*  0  0 


(1.8) 


(1.9) 


(1.10) 


By  Ampere's  law,  the  azimuthal  magnetic  field  is  related  to  the  anode 
current  by  the  expression, 


(1.11  ) 
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whereas  the  radial  electric  field  is  given  by  Equation  (1.5) 
of  Equations  (1.5)  and  (1.11)  in  (1.8),  we  find  that 


—  =  <ue) 

a  0  O 

B9 


-1/2 


Making  use 


1.12) 


as  expected. 


V>>1 


Figure  1.2  presents  simulation  results  for  the  inlet  voltage  as 
a  function  of  time.  Note  that  the  imposition  of  the  boundary  condition  <t>” 
=  0  forces  the  total  wave  to  be  exactly  equal  to  the  incident  wave,  or  $  = 
4>+.  The  voltage  at  the  outlet  is  presented  in  Figure  1.3.  Note  the 
following:  (1)  the  time  delay  correctly  represents  a  wave  velocity  c,  (2) 
the  outlet  voltage  is  invariant  ( 4>  =  10 6  V)  at  late  time,  and  (3) 
formation  of  a  foot  near  the  leading  edge  and  a  ripple  near  the  peak. 
This  result,  also  in  evidence  for  the  outlet  magnetic  field  in  Figure  1.4, 
is  a  consequence  of  numerical  diffusion  in  the  code. 

Figures  1.5,  1.6,  and  1.7  present  electric  and  magnetic  field 
distributions.  All  illustrate  the  r” 1  behavior  of  the  TEM  wave  predicted 
in  Equation  (1.5).  The  result  in  Figure  1.5  was  imposed  as  a  boundary 
condition;  however,  the  very  similar  result  for  the  outlet  field  in  Figure 
1.6  is  not  so  constrained.  Comparison  suggests  the  potential  accuracy  of 
the  method  in  purely  electromagnetic  applications. 


As  a  final  check,  we  compute  using  Equations  (1.3)  and  (1.11) 
the  magnetic  fields  at  r  =  0.105  m  and  0.195  m.  (Note  that  the  azimuthal 

magnetic  field  is  defined  at  half-grid  rather  than  full-grid  points). 

2  2 

The  resulting  analytical  results,  B  =  4.58x10”  and  2.47x10”  T,  are  given 
exactly  (to  three  significant  figures)  by  the  code. 
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Table  1.1.  Input  data  for  Problem  1.1. 


title  *  problem  1.1  *  / 
comment  *  base  problem  *  / 
system  2  / 

xlgrid  1  62  2  0.0  60  0.01  0.6  / 
x2grid  1  12  2  0.1  10  0.01  0.1  / 
fields  1  1  500  2.0&-11  / 
courant  00/ 

conductor  cathode  12  2  62  2  / 

conductor  anode  -1  2  12  62  12  / 

voltage  temporal  radial  twod  000222  12  2/ 

function  temporal  0  3  0.0  0.0  1.0e-9  1.0e6  1.0  1.0e6  / 

function  radial  5  -11/ 


lookback  1  -1  62  2  62  12  / 
observe  1  twod  10.01.02  22  2  12 

twod  1  0.0  1.0  2  32  2  32  12 

twod  1  0.0  1.0  2  62  2  62  12 

twod  10. 01. 06  22  22 

twod  10. 01. 06  32  2  32  2 
twod  10. 01. 06  61  2  61  2 
twod  10. 01. 062  11  2  11 
twod  1  0.0  1.0  6  32  11  32  11 

twod  1  0.0  1.0  6  61  11  61  11  / 


range  250  12  22  2  12  1 

1  2  32  2  32  12  1 

1  2  62  2  62  12  1 

1  6  2  2  2  12  1 

1  6  32  2  32  12  1 

1  6  61  2  61  12  1  / 


output  0  / 
start  / 
stop  / 
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MAGIC  VERSION:  JUNE  1933  GFTE :  34/07/02 
SIMULATION:  PROBLEM  1.1 


TIME  HISTORY  PLOT 
E2  COMPONENT 

INTEGRATED  FROM  (2,2)  TO  (2,12) 


4.0 

TIME 


Figure  1.2.  Inlet  voltage  vs.  time,  Problem  1.1. 
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1.2  SCATTERED  WAVE  REFLECTION  AT  INLET. 

1.2.1  Problem  Description. 

This  problem  demonstrates  the  application  of  the  previously 
described  (zero  scattered  wave)  inlet  boundary  condition  for  a  case  in  in 
which  there  is  a  physically  meaningful  scattered  wave  impinging  on  it, 
i.e.,  t>~  *  0.  To  create  a  scattered  wave,  we  shall  simply  short  out  the 
two  conductors  at  the  previous  outlet.  This  should  produce  a  scattered 
wave  equal  in  magnitude  to  the  Incident  wave  ($'  =  -4>+).  In  this  problem, 
we  will  see  that  the  selected  inlet  boundary  condition  ( =  0)  is  the 
wrong  boundary  condition  for  this  case. 

1.2.2  Suggested  Approach. 

Using  the  input  data  for  Problem  1.1,  simply  replace  the 
LOOKBACK  command  at  the  outlet  with  a  CONDUCTOR  command.  This  is  labeled 
'SHORT'  in  TabLe  1.2. 

1.2.3  Analytical  Solution. 

To  understand  (or  predict)  what  will  happen  in  this  problem,  it 
is  necessary  to  understand  the  effect  of  an  impedance  mismatch.  Thus,  we 
consider  a  two-region  problem,  with  the  inlet  at  impedance  Z^  and  the 
outlet  at  Z2. 

As  before,  the  incident  pulse  is  defined  by  a  peak  voltage,  i>+, 
and  associated  current,  I+,  which  are  related  by  the  inlet  impedance 
according  to 

=  I-  Zl  .  (1.13) 


1-16 


At  the  mismatch,  scattering  occurs.  The  scattered  wave,  described  in 
terms  of  peak  voltage  and  current,  and  I',  propagates  in  both 

directions  from  the  discontinuity  at  the  speed  of  light.  The  voltage- 
current  relationship  must  be  satisfied  on  both  sides  of  the  discon¬ 
tinuity.  Thus,  solving  for  the  scattered  voltage  and  current  in  terms  of 
the  incident  pulse,  one  obtains 


z2  - 

Zl 

Zi  + 

Z2 

.  Zl  - 

12 

Zi  ♦ 

ll 

iltage 

and 

2Z2 

Zi  ♦  i 

-2 

2Zi 

(1.14; 


Then  the  total  voltage  and  current  are  given  by 


(1.15) 


Zi  ♦  Z2 


For  a  short  circuit  (Z2  =  0).  the  total  voltage  vanishes  and  the  current 
doubles.  Note  that  the  transient  results  should  reflect  the  retarded  time 
between  mismatch  and  measurements  as  well  as  the  rise-time  of  the  incident 
pulse. 


Figure  1.9  illustrates  the  simulation  result  at  the  inlet 
boundary.  Note  that  imposition  of  the  condition  $'  =  0  forces  the  total 
wave  to  be  equal  to  the  Incident  wave,  or  $  =  $+.  Physically,  we  expect 
the  voltage  at  the  inlet  to  vanish  at  5  nsec. 
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The  reason  for  the  discrepancy  is  clear  in  Figure  1.9,  which 
illustrates  voltage  vs.  time  at  midpoint  (L/2).  Note  the  passage  of  the 
incident  wave  beginning  at  1  nsec.  At  3  nsec,  the  scattered  wave  from 
the  short  circuit  (traveling  from  right  to  left)  begins  to  cancel  the 
incident  wave.  Between  4  and  5  nsec,  cancellation  is  complete 
( <t>+  +  4>"  =  0)  and  the  voltage  vanishes.  So  far,  the  result  is  physically 
correct.  However,  by  5  nsec  the  scattered  wave  from  the  short  circuit 
has  reached  the  inlet,  has  itself  been  rescattered  (since  =  O', 

and  the  second  scattered  wave  has  returned  to  the  midpoint.  This  process 
will  simply  repeat  forever  -  of  course,  we  have  picked  the  spatial 
dimensions  and  wave  rise-time  to  produce  the  even  pattern  shown  in  Figure 
1 .9. 

This  result  is,  of  course,  nonphysical  and  potentially 
catastrophic.  This  is  illustrated  even  more  clearly  in  the  plot  of 

magnetic  field  vs.  time.  This  result,  shown  in  Figure  1.10,  indicates 

that  the  magnetic  field  will  Increase  without  limit.  (This  result  is 

similarly  predicted,  with  careful  attention  to  signs.)  Thus,  it  is 
crucially  important  for  simulation  purposes  to  separate  experimental 
measurements  into  incident  and  scattered  waves.  Application  of  the  total 
wave  will  produce  satisfactory  results  only  if  the  choice  of  boundary 
condition  is  consistent  with  the  internal  dynamics.  Otherwise,  scattered 
waves  will  be  trapped  in  the  simulation.  In  Problem  1.3,  we  consider  a 
boundary  condition  which  avoids  this  particular  difficulty. 


Table  1.2. 


Input  data  for  Problem  1.2. 


title  *  problem  1.2  *  / 

comment  *  scattered  wave  reflection  *  / 
system  2  / 

xlgrid  1  62  2  0.0  60  0.01  0.6  / 

x2grid  1  12  2  0.1  10  0.01  0.1  / 

fields  1  1  500  2.0e-ll  / 
courant  00/ 

conductor  cathode  12  2  62  2  / 
conductor  anode  -1  2  12  62  12  / 
conductor  short  -1  62  2  62  12  / 

voltage  temporal  radial  twod  000222  12  2/ 
function  temporal  0  3  0.0  0.0  1.0e-9  1.0e6  1.0  1 . 0e6  / 
function  radial  5  -11/ 

observe  1  twod  10.01.02  22  212 

twod  1  0.0  1.0  2  32  2  32  12 

twod  1  0.0  1.0  2  62  2  62  12 

twod  10. 01. 06  22  22 

twod  1  0.0  1.0  6  32  2  32  2 

twod  1  0.0  1.0  6  61  2  61  2 

twod  10. 01. 062  11  2  11 
twod  1  0.0  1.0  6  32  11  32  11 


twod  1  0.0  1.0  6  61  11  61  11  / 
range  250  12  22  2  12  1 

1  2  32  2  32  12  1 

1  2  62  2  62  12  1 

1  6  2  2  2  12  1 

1  6  32  2  32  12  1 

1  6  61  2  61  12  1  / 

output  0  / 
start  / 
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1.3  SCATTERED  WAVE  ABSORPTION  AT  INLET. 

1.3.1  Problem  Description. 

The  problem  is  to  repeat  the  short  circuit  simulation  in  Problem 
1.2,  but  to  use  an  inlet  boundary  condition  which  allows  the  scattered 
waves  to  exit  through  the  iniet  rather  than  forcing  them  to  zero.  This  is 
accomplished  by  absorbing  the  scattered  wave  at  the  inlet. 

1.3.2  Suggested  Approach. 

The  input  data  from  Problem  1.2  can  be  duplicated  exactly  except 
for  the  VOLTAGE  command.  Here,  the  option  flag  (which  doubles  as  a 

surface  alignment  flag)  should  be  set  to  +1;  also  the  phase  velocity 
factor  should  be  set  to  unity.  The  input  data  is  shown  in  Table  1.3. 

1.3.3  Analytical  Solution. 

Selected  results  from  this  simulation  are  shown  in  Figures  1.11. 
1.12,  and  1.13,  which  correspond  to  similar  plots  from  Problem  1.2. 

Whereas,  in  Figure  1  .8  the  total  voltage  is  maintained  at  the  incident 

wave  value  (i  =  i>+),  In  Figure  1.11  the  scattered  wave  is  correctly 

accounted  for  ( $  =  $+  +  $*).  Note  that  the  time  lapse  until  the  voltage 
begins  to  drop  is  perfectly  consistent  with  the  leading  edge  of  the 

incident  pulse  reaching  the  short  circuit,  scattering,  and  returning  to 
the  inlet. 
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The  voltage  and  current  near  midpoint  (L/2)  are  shown  in  Figures 
1.12  and  1.13.  From  Equation  (1.15),  we  expect  the  voltage  to  vanish  and 
the  current  to  double,  precisely  as  the  figures  illustrate.  3gain,  note 
the  consistency  of  the  time  lapses  with  waves  traveling  at  the  speed  of 
light . 
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In  summary,  this  inlet  boundary  condition,  which  allows 
scattered  waves  to  escape  upstream  (absorption),  is  generally  superior  to 
that  used  in  Problems  1.1  and  1.2.  However,  it  suffers  from  two  defects: 
(1)  the  phase  velocity  for  the  scattered  wave  must  be  known  (see  Problem 
1.5),  and  (2)  the  incident  wave  must  be  clearly  distinguishable  from  any 
scattered  waves. 
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SCATTERED  WAVE  ABSORPTION  AT  INLET  AND  OUTLET. 


1  .4 

1.4.1  Problem  Description. 

Repeat  the  previous  simulation,  allowing  scattered  waves  to  be 
absorbed  at  both  inlet  and  outlet. 

1.4.2  Suggested  approach. 

The  input  data  from  Problem  1.3  can  be  duplicated  exactly, 
except  for  the  CONDUCTOR  command  which  specifies  the  short.  This  should 
be  replaced  with  a  LOOKBACK  command,  as  shown  in  Table  1.4. 

1.4.3  Analytical  Solution. 

Results  from  this  simulation,  as  shown  in  Figures  1.14  and  1.15, 
closely  duplicate  those  from  Problem  1.1.  Note,  however,  that  the  inlet 
voltage  in  Figure  1.14  does  not  exactly  reproduce  the  incident  wave,  due 
to  a  very  small  but  non-zero  scattered  component  originating  from  numeri¬ 
cal  diffusion. 


Table  1.4.  Input  data  for  Problem  1.4. 


title  *  problem  1.4  *  / 

comment  *  lookback  reinstated  *  / 

system  2  / 

xlgrid  1  62  2  0.0  60  0.01  0.6  / 
x2grid  1  12  2  0.1  10  0.01  0.1  / 
fields  1  1  500  2.0e-ll  / 
courant  00/ 

conductor  cathode  12  2  62  2  / 
conductor  anode  -1  2  12  62  12  / 
voltage  temporal  radial  twod  01. 01222  12  2 
function  temporal  0  3  0.0  0.0  1.0e-9  1 . 0e6  1.0 


function 

radial  5 

1 

-11/ 

lookback 

1  ■ 

-1  62 

2 

62  12  / 

observe 

1 

twod 

1 

0.0  1.0 

2 

2 

2 

2 

12 

twod 

1 

0.0  1.0 

2 

32 

2 

32 

12 

twod 

1 

0.0  1.0 

2 

61 

2 

61 

12 

twod 

1 

0.0  1.0 

6 

2 

2 

2 

2 

twod 

1 

0.0  1.0 

6 

32 

2 

32 

2 

twod 

1 

0.0  1.0 

6 

61 

2 

61 

2 

twod 

1 

0.0  1.0 

6 

2  : 

11 

2  11 

twod 

1 

0.0  1.0 

6 

32 

11 

32 

11 

twod 

1 

0.0  1.0 

6 

61 

11 

61 

11  / 

range  250 

1 

2  2 

2 

2  12  1 

1 

2  32 

2 

32  12  1 

i_ 

2  61 

2 

61  12  1 

1 

6  2 

2 

2  12  1 

1 

6  32 

2 

32  12  1 

1 

6  61 

2 

61  12  1 

/ 

/ 

1 . 0e6  / 


output  0  / 
start  / 


Table  1.5.  Input  data  for  Problem  1.5. 


title  *  problem  1.5  *  / 

comment  *  lookback  phase  velocity  *  / 

system  2  / 

xlgrid  1  62  2  0.0  60  0.01  0.6  / 
x2grid  1  12  2  0.1  10  0.01  0.1  / 
fields  1  1  500  2.0e-ll  / 
courant  00/ 

conductor  cathode  12  2  62  2  / 
conductor  anode  -1  2  12  62  12  / 

voltage  temporal  radial  twod  01. 01222  12  2/ 
function  temporal  0  3  0.0  0.0  1.0e-9  1.0e6  1.0  1.0e6  / 
function  radial  5  -11/ 

lookback  0.8  -1  62  2  62  12  / 
observe  1  twod  10. 01. 02  22  2  12 

twod  1  0.0  1.0  2  32  2  32  12 

twod  1  0.0  1.0  2  61  2  61  12 

twod  10. 01. 06  22  22 

twod  10. 01. 06  32  2  32  2 
twod  1  0.0  1.0  6  61  2  61  2 
twod  10. 01. 062  11  2  11 
twod  1  0.0  1.0  6  32  11  32  11 

twod  1  0.0  1.0  6  61  11  61  11  / 

range  250  12  22  2  12  1 

1  2  32  2  32  12  1 

1  2  61  2  61  12  1 

1  6  2  2  2  12  1 

1  6  32  2  32  12  1 

1  6  61  2  61  12  1  / 

output  0  / 
start  / 


vs 


COAXIAL  LINE  WITH  MISMATCH. 


1  .6. 1 


Problem  Description. 


The  preceding  problems  have  explored  the  effects  of  various 
boundary  conditions  and  artificial  impedance  effects  which  can  result  in 
simulations.  In  this  problem,  we  consider  a  physical  mismatch:  the  anode 
radius  is  to  be  reduced  from  0.20  m  to  0.15  m  from  the  axial  midpoint 
(L/2)  onward.  This  coaxial  line,  illustrated  in  Figure  1.18,  has  the 
dimensions, 


=  0.10  m 
=  0.15  m 
=  0.20  m 
=  0.60  m 


(1.17) 


1.6.2  Suggested  Approach. 

We  suggest  using  the  spatial  grid,  etc.,  from  the  preceding 
simulations,  with  a  few  of  the  boundaries  and  measurements  altered.  For 
example,  the  CONDUCTOR  command  for  the  anode  surface  can  be  altered  to 
produce  the  desired  geometry  using  the  continuation  option.  (The  input 
data  for  the  first  line  segment  is  produced  in  its  entirety,  including 
alphanumeric,  surface  normal  flag,  and  spatial  line  indices.  Only  the  end 
point  indices  are  specified  for  the  remaining  contiguous  line  segments.) 
In  addition,  the  spatial  line  indices  for  the  LOOKBACK  command  and  various 
output  commands  must  be  modified  to  reflect  the  smaller  radius.  The  input 
data  for  this  problem  is  shown  in  Table  1.6. 


1  .6.3 


Analvtical  Solution. 


This  is  a  classic  example  of  the  two-region  impedance  mismatch. 
From  Equations  (1.9)  and  (1.17),  the  impedances  of  the  inlet  and  outlet 
regions  are 


V  *.*  s.'  •;  •.'.v.x  .s  .v  v  v  V  V  V  V  V  V  v  •.*. 
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Zi  =  41.56  a 

( 1 . IS) 

Z2  =24.31  a  , 

respectively.  Making  use  of  Equation  (1.15),  we  obtain  the  total 
steady-state  voltage  and  current, 

D  =  7 . 39x10  5  V 

4  (1.19) 

I  =  3.03x10  A  . 


The  simulation  results  for  inlet  and  outlet  voltage  are  shown  in 
Figures  1.19  and  1.20,  respectively,  with  the  corresponding  magnetic 
fields  in  Figures  1.21  and  1.22.  By  applying  Equation  (1.11)  to  the 
magnetic  field  (both  voltage  and  field  data  time  averaged  over  final  1 
nsec),  we  obtain  the  simulation  results, 


<D  =  7.38x10 5  V 
I  =  3.04x104  A  , 


(1.20) 


which  compare  favorably  with  the  analytical  results.  The  electric  and 
magnetic  field  distributions  (vs.  radius)  should  also  be  checked  to 
confirm  the  existence  of  TEM  waves  at  inlet  and  outlet. 
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Table  1.6.  Input  data  for  Problem  1.6. 


title  *  problem  1.6*/ 
comment  *  impedance  mismatch  *  / 
system  2  / 

xlgrid  1  62  2  0.0  60  0.01  0.6  / 
x2grid  1  12  2  0.1  10  0.01  0.1  / 
fields  1  1  500  2.0e-ll  / 
courant  00/ 

conductor  cathode  12  2  62  2  / 
conductor  anode  -1  2  12  32  12  32  7  62  7  / 
voltage  temporal  radial  twod  01. 01222  12  2 
function  temporal  0  3  0.0  0.0  1.0e-9  1.0e6  1.0 


function 

radial  5 

-11/ 

lookback 

1  - 

-1  62 

2 

62  12  / 

observe 

1 

twod 

1 

0.0  1.0 

2 

2  2 

2 

12 

twod 

1 

0.0  1.0 

2 

31  2 

31 

12 

twod 

1 

0.0  1.0 

2 

61  2 

61 

7 

twod 

1 

0.0  1.0 

6 

2  2 

2 

2 

twod 

1 

0.0  1.0 

6 

31  2 

31 

2 

twod 

1 

0.0  1.0 

6 

61  2 

61 

2 

twod 

1 

0.0  1.0 

6 

2  11 

2 

11 

twod 

1 

0.0  1.0 

6 

31  11 

.  31  11 

twod 

1 

0.0  1.0 

6 

61  6 

61 

6  / 

range  250 

1 

2  2 

2 

2  12  1 

1 

2  31 

2 

31  12  1 

1 

2  61 

2 

61  7  1 

1 

6  2 

2 

2  12  1 

1 

6  31 

2 

31  12  1 

1 

6  61 

2 

61  7  1 

/ 

output  0  / 
start  / 
stop  / 


/ 

1 . 0e6  / 


1  _"50 


Axis  of  Symmetry 


on  01 


coaxial 


1  i ne  wi th 


mi smatcn . 
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1 .7 


REDUCED  TIME  STEP. 


M 


1.7.1  Problem  Description. 

We  want  to  explore  the  electromagnetic  effects  of  using  a 
smaller  time  step,  well  below  the  Courant  criterion. 

1.7.2  Suggested  Approach. 

Using  input  data  from  the  preceding  simulation,  cut  the  time 
step  in  the  FIELDS  command  in  half.  Also,  double  the  number  of  time  steps 
and  all  time  step  intervals  that  determine  the  frequency  of  output.  The 
new  input  data  is  shown  in  Table  1.7. 


Vjv  1 
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1.7.3  Analytical  Solution. 

Figure  1.23  illustrates  the  calculated  inlet  magnetic  field  vs. 
time.  This  is  one  of  the  more  sensitive  results,  and  can  be  compared  with 
the  result  of  Figure  1.21.  While  there  are  slight  discernible 
differences,  the  agreement  is  substantial.  This  is  a  characteristic  of 
the  centered-dif ference  electromagnetic  algorithm  -  results  generally 
depend  upon  spatial  differences  but  not  upon  temporal  differences. 


Table  1.7.  Input  data  for  Problem  1.7. 


title  *  problem  1.7  *  / 
comment  *  reduced  time  step  *  / 
system  2  / 

xlgrid  1  62  2  0.0  60  0.01  0.6  / 
x2grid  1  12  2  0.1  10  0.01  0.1  / 
fields  1  1  1000  1.0e-ll  / 
courant  00/ 

conductor  cathode  12  2  62  2  / 
conductor  anode  -1  2  12  32  12  32  7  62  7  / 
voltage  temporal  radial  twod  01. 01222  12  2/ 
function  temporal  0  3  0.0  0.0  1.0e-9  1.0e6  1.0  1.0e6  / 


function 

radial  5 

l 

-11/ 

lookback 

1  - 

-1  62 

2 

62  12  / 

observe 

1 

twod 

1 

0.0  1.0 

2 

2  2 

2 

12 

twod 

1 

0.0  1.0 

2 

31  2 

31 

12 

twod 

1 

0.0  1.0 

2 

61  2 

61 

7 

twod 

1 

0.0  1.0 

6 

2  2 

2 

2 

twod 

1 

0.0  1.0 

6 

31  2 

31 

2 

twod 

1 

0.0  1.0 

6 

61  2 

61 

2 

twod 

1 

0.0  1.0 

6 

2  11 

2 

11 

twod 

1 

0.0  1.0 

6 

31  11  31  11 

twod 

1 

0.0  1.0 

6 

61  6 

61 

6  / 

range  250 

1 

2  2 

2 

2  12  1 

1 

2  31 

2 

31  12  1 

1 

2  61 

2 

61  7  1 

1 

6  2 

2 

2  12  1 

1 

6  31 

2 

31  12  1 

1 

6  61 

2 

61  7  1 

/ 

output  0  / 
start  / 
stop  / 


1-46 


TIME-BIASED  ALGORITHM. 


1.9.1  Problem  Description. 

As  a  final  exercise,  we  wish  to  duplicate  Problem  1.6  using  the 
time-biased  electromagnetic  algorithm.  This  algorithm,  which  is  implicit, 
is  most  useful  in  damping  out  the  numerical  noise  associated  with  particle 
kinematics. 

1.8.2  Suggested  Approach. 

From  the  User's  Manual,  we  arbitrarily  select  time-biased 
coefficients  for  four  iterations, 


ai  =  a2  =  0.5  ,  a3  =  0 
t  =  0.29912,  0.15022,  0.11111,  U1,. 


(1.21) 


From  the  Courant  criterion, 


1  2  ,  1  1  , 

- - -  >  (a2  -  4a1a3)  {  - -  ♦  - -  }  , 

:26t 2  ( 6r)  2  (6z)2 


(1.22) 


it  is  obviously  possible  to  double  the  time  step  and  to  cut  the  number  of 
time  steps  and  output  frequency  in  half.  (Note  that  the  running  time  will 
actually  Increase  by  about  100%;  the  time  step  is  doubled,  but  four 
iterations  are  specified  per  step.)  The  input  data  for  this  simulation  is 
given  in  Table  1.8. 

1.8.3  Analytical  Solution. 

The  calculated  magnetic  field  at  the  inlet  is  shown  in  Figure 
1.24.  This  result  can  be  compared  with  Figures  1.21  and  1.23.  While  the 
mean  results  are  in  good  agreement,  the  result  in  Figure  1.24  is 
relatively  free  of  high-frequency  (dispersion)  numerical  noise.  This 
ability  to  damp  high-frequency  waves  Is  a  primary  characteristic  of  the 
time-biased  algorithm. 
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Table  I.S.  Input  data  for  Problem  1.5. 


title  *  problem  1.8  *  / 

comment  *  time-biased  field  algorithm  *  / 
system  2  / 

xlgrid  1  62  2  0.0  60  0.01  0.6  / 
x2grid  1  12  2  0.1  10  0.01  0.1  / 
fields  1  3  250  4.0e-ll  0.5  0.5  0.0 

4  1.0  0.29912  0.15022  0.11111  / 

courant  00/ 

conductor  cathode  12  2  62  2  / 
conductor  anode  -1  2  12  32  12  32  1  62  1  / 
voltage  temporal  radial  twod  01.01222122/ 
function  temporal  0  3  0.0  0.0  1.0e-9  1.0e6  1.0  1.0e6  / 


function  radial  5 

» 

-11/ 

lookback  1  ■ 

-1  62 

2 

62  12  / 

observe  1 

twod 

1 

0.0  1.0 

2 

2  2 

2 

12 

twod 

1 

0.0  1.0 

2 

31  2 

31 

12 

twod 

1 

0.0  1.0 

2 

61  2 

61 

7 

twod 

1 

0.0  1.0 

6 

2  2 

2 

2 

twod 

1 

0.0  1.0 

6 

31  2 

31 

2 

twod 

1 

0.0  1.0 

6 

61  2 

61 

2 

twod 

1 

0.0  1.0 

6 

2  11 

2 

11 

twod 

1 

0.0  1.0 

6 

31  11 

31  11 

twod 

1 

0.0  1.0 

6 

61  6 

61 

6  / 

range  125  1 

2  2 

2 

2  12  1 

1 

2  31 

2 

31  12  1 

1 

2  61 

2 

61  7  1 

1 

6  2 

2 

2  12  1 

1 

6  31 

2 

31  12  1 

1 

6  61 

2 

61  7  1 

/ 

output  0  / 
start  / 
stop  / 
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SECTION  2 

PARTICLE  KINEMATICS 

The  problem  set  on  particle  kinematics  is  introduced  by  simula¬ 
tions  involving  the  motion  of  a  single  particle  in  a  prescribed  electro¬ 
magnetic  field.  Four  problems  were  chosen  to  represent  various  physical 
effects,  including  relativistic  motion,  gyromagnetic  motion,  E  x  B  drift, 
and  3  x  VB  drift.  These  single  particle  problems  also  demonstrate  the 
extent  to  which  the  code  can  faithfully  reproduce  physical  particle 
trajectories.  Following  the  single-particle  problems  in  this  section  are 
two  fully  self-consistent  problems  in  one-dimension.  These  two  problems 
cover  field  emission  and  magnetic  insulation. 
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To  test  the  particle  kinematics  algorithm,  KAGIC  allows  the  user 
to  prescribe  external  fields  causing  forces  on  the  particles.  The  FORCES 
command  allows  the  prescription  of  any  (or  all  six)  field  components  as  a 
triple  product  function  in  space  and  time.  For  problem  2.1,  all  fields 
except  E2  are  zero,  and 

E2(xi,X2,t)  =  f  1  (x  i )  f  2(x  2)  ft(t)  .  (2.5) 

The  functions  are  identified  by  means  of  an  option  in  the  FORCES  data 
sequence.  Each  separate  function  is  specified  in  the  usual  way  with  a 
FUNCTION  command.  This  capability  provides  a  simple  means  of  prescribing 
fields  required  for  a  particular  kinematics  test,  e.g.,  8  x  7B,  which 
clearly  requires  a  spatial  dependence. 

When  prescribed  fields  are  used,  it  is  still  necessary  to 
provide  input  data  for  all  the  usual  data  sequences,  including  the  field 
algorithm  and  spatial  grid.  We  suggest  using  a  1  x  1  x  1  mesh  (i.e.,  one 
functional  cell  sufficiently  large  to  contain  the  particle  trajectory) 
surrounded  by  conducting  boundaries.  In  this  case,  the  field  algorithm  is 
essentially  ignored,  and  the  time  step  specified  by  the  FIELDS  command  is 
used  only  for  the  particle  kinematics.  Output  for  these  simulations  can 
easily  be  obtained  in  the  form  of  trajectory  plots,  which  can  be  measured 
directly  to  verify  analytical  results.  More  precise  measurements  are 
available  from  the  diagnostics  which  print  out  particle  coordinates  and 
momenta  during  the  simulation. 

This  problem  is  easily  set  up  by  including  three  commands: 
FORCES  with  the  prescribed  field  option  to  set  the  electric  field. 
POPULATE  to  establish  Initial  conditions  for  a  single  particle,  and 
DIAGNOSE  (KINEMAT)  to  print  particle  coordinates  as  a  function  of  time. 


(A  trajectory  plot  should  be  specified  with  the  TRAJECTORY  command;  how¬ 
ever,  it  will  not  allow  quantitative  measurements.  The  object  is  to 
obtain  x  and  y  as  a  function  of  t,  and  to  compare  computed  with  analytical 
results  at  several  points  in  time.)  The  input  data  for  this  simulation  is 
shown  in  Table  2.1. 

2.1.3  Analytical  Solution. 

This  problem,  although  relativistic,  is  analytical  (inte- 
grable).  The  Lagrangian  is  given  by^ 


L  =  -  mc2(l  - 


•2  • 2  1/2 
x  +  y  'i  ' 


)  -  eEy 


Application  of  the  variational  principle, 


d  fdL-i  5L 

^  V  '  ^ 


=  0  , 


(2.6) 


(2.7) 


and  integration  yield  the  results. 


•2  "2  1/2 

(i  - 


(2.8) 


•2  *2  1/2 
x  *  y  ) 

O  * 


*  t  . 
m 


By  algebraic  manipulation,  we  obtain  the  two  integrals 


=  P  i 


„  2  -  ,  2  1/  2 
{[,.(£)  ]  .  (-  ±  t  )  | 


t  H.  Goldstein,  Classical  Mechanics,  Addison-Wesley ,  New  York,  1950 

(p.206) . 
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(2.9) 


y  = 


(-  —  t  )  dt 
m 


{  [1  ♦  (P)  ]  ♦  (-  *  t 

c  J  v  me 


eE  '^2M2 


These  are  integrated  to  give  the  trajectory  equations, 


-  ±  t  ♦  [i  +  (?)  +  (£1  t) 

me  L  c  me 


2  eE 


X  = 


mo  «  f 

—  P  In  1 
eE 


[1  -  <?>  ] 


2  1/2 


(2.10) 


y  = 


me 


eE 


2  an  2  l/2 

-  eE  ^  \  i  < 


{  [1  -  (&  +  (||  t)  ]  -  [1  ♦  (?)  ] 


2  V2, 


Although  there  is  no  force  in  the  x  direction,  this  component  of 
velocity,  since  initially  non-zero,  must  decrease  due  to  the  increase  in 
mass.  Thus,  x  is  sublinear  in  time.  This  effect  appears  in  the 
simulation  trajectory  plot  illustrated  in  Figure  2.1.  Examining 
diagnostic  results  from  the  output  listing,  we  find  the  calculated 
coordinates  at  t  =  0.141  nsec  to  be 


x  =  8.97x10"3  m 


y  =  2.94x10"  m 


(2.11  ) 


The  analytical  result  from  Equation  (2.10)  at  the  same  time  yields 


x  =  9.06x10"3  m 
y  =  2 . 81x10" 2  m  . 


(2.12) 


This  time  period  includes  only  50  time  steps  in  the  kinematics:  the 
accuracy  of  the  numerical  result  will  improve  with  resolution. 
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Table  2.1.  Input  data  for  Problem  2.1. 


title  *problem  2.1*  / 

comment  ^relativistic  trajectory  in  constant  electric  field*  / 
system  1  / 

xlgrid  132  0.0  1  0.028  0.028  / 

x2grid  132  0.0  1  0.028  0.028  / 

conductor  boundary  +12232332322/ 

fields  1  1  50  2.82616e-12  / 

courant  00/ 

kinematics  10000/ 

forces  1  1  1  2  e2  null  null  / 

function  e2  1  0  -3.0e7  / 

populate  electrons  00. 012233  1.6e-19 

0.0  0.0  1.0e8  0.0  0.0  / 

diagnose  prscrib  0  2 

kinemat  0  51  / 

trajectory  50  50  1  0.0  0.028  0.0  0.028  / 
timeout  -10/ 
start  / 
stop  / 
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GYROMAGNETIC  MOTION. 


2.2.1  Problem  Description. 


This  problem  involves  simulating  the  trajectory  of  a  relativis¬ 
tic  electron  in  a  high  magnetic  field, 


B  =  zB  , 


where  B  =  -  20  T.  Initial  conditions  for  the  particle  are 


(2.13) 


t  =  x  =  y  =  z  =  0 


P  =  xp 


(2.14) 


where  p  =  9  x  10B  m/sec. 


2.2.2  Suggested  Approach. 


In  this  problem,  the  trajectory  plot  capability  is  appropriate 
to  measure  performance.  One  trick  is  to  choose  the  kinematics  time  step 
such  that  an  integer  multiple  equals  the  gyromagnetic  period,  i.e., 


k5t  =  t  . 


(2.15) 


Then  the  orbit  should  close  in  k  time  steps.  Conservation  of  momentum, 
etc.,  can  be  assessed  using  the  previously  described  diagnostic 
capability.  Input  data  for  this  simulation  is  shown  in  Table  2.2. 

2.2.3  Analytical  Solution. 

In  the  absence  of  electric  field  (E  =  0),  the  momentum  magnitude 
(and  thus  y)  is  constant.  The  gyromagnetic  radius  and  period  are  given  by 
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For  the  momentum  specified  in  Equation  (2.14),  the  radius  and  period  are 
0.2558  mm  and  5.65  psec,  respectively. 

It  will  be  seen  from  the  simulation  trajectory  plot  in  Figure 
2.2  that  the  orbit  has  apparently  closed  in  100  time  steps.  Using  the 
kinematics  diagnostic  data  from  the  output  listing,  we  estimate  the 
effective  radius  to  be  0.2559  mm.  This  degree  of  accuracy  results  from 
using  a  small  time  step  (i.e.,  large  number  of  steps  to  represent  a  single 
orbit).  With  a  larger  time  step,  the  accuracy  decreases  and,  in  the 
limit,  approaches  a  simple  oscillating  state. 


>  ", J>  1*  w  n* 


■j.  vv  -av:#  -jtj>  vvvA'A>v.vrATAV.v;/r/rA'AV.|,-V.,;/.v..- 


Table  2.2.  Input  data  for  Problem  2.2. 


title  *problem  2.2*  / 

comment  *gyromagnetic  motion  in  constant  magnetic  field*  / 
system  1  / 

xlgrid  132  -5.0000e-4  1  0.001  0.001  / 

x2grid  132  -7.5587e-4  1  0.001  0.001  / 

conductor  boundary  +12232332322/ 

fields  1  1  100  5 . 65233e-14  / 

courant  00/ 

kinematics  10000/ 

forces  1116  b3  null  null  / 

function  b3  1  0  -20.0  / 

populate  electrons  00. 012233  1.6e-19 


0.5  0.75587  9.0e8  0.0  0.0  / 
diagnose  prscrib  0  2 

kinemat  0  101  / 

trajectory  100  100  1  -0.0005  0.0005  — 7.5587e— 4  2.44l3e-4  / 
timeout  -10/ 
start  / 
stop  / 
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Figure  2.2.  Trajectory  plot,  Problem  2.2 
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2.3 


DRIFT  (E  x  B). 


2.3.1 


Problem  Description. 


This  problem  involves  simulating  an  E  x  B  drift  where  the 
constant  fields  obey  the  constraint, 


—  <  1  . 

'cB 


(2.17) 


This  constraint  is  motivated  by  the  reduction  in  normal  electric  field 
which  results  from  space-charge  limiting  in  magnetic  insulation.  We 
arbitrarily  select  field  parameters  from  the  preceding  problems,  i.e., 

®  =  ZB  (2.18) 
E  =  yE  , 

where  B  =  -20  T  and  E  =  -3  x  10 7  V/m.  The  particle  initial  conditions  are 


2.3.2 


t  =  x  =  p  =  0  . 


Suggested  Approach. 


(2.19) 


Using  a  trajectory  plot  (and  diagnostic),  measure  the  period, 
turning  distance,  and  drift  of  the  orbit.  Input  data  for  the  simulation 
is  shown  in  Table  2.3. 


2.3.3 


Analytical  Solution. 


For  the  specified  fields  (which  satisfy  the  constraint  of 
Equation  (2.17)),  a  Lorentz  transformation  can  be  made  to  a  frame,  s', 
in  which  the  electric  field  vanishes.  This  frame,  illustrated  in  Figure 
2.3,  is  defined  by  the  velocity, 


E 

B 


(2.20) 
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The  only  nonvanishing  field  in  the  5  frame  is  the  magnetic  field,  which 
transforms  according  to 


B  = 

.  1 


Y~  XB 


2  1/2 

(1  - 


[1  -  £-)  ] 

cB 


2  l/2 


(2.21  ) 


As  illustrated  in  Figure  2.3b,  the  turning  distance  in  this  frame  is 
related  to  gyromagnetic  radius  by 


d'  =  2r'  = 


2p'  _  2m£y" 


(2.22) 


eB'  eB‘ 


The  period  and  drift  are  given  by 

2 


2ur'  2rany 


eB 


(2.23) 


V  =  0  . 


These  parameters  must  be  transformed  back  to  the  rest  frame.  Since  the 

perpendicular  dimension  is  Lorentz  invariant,  the  turning  distance  is  the 

same  in  either  frame  (i.e.,  d  =  d').  However,  the  period  and  drift  are 

modified.  Thus,  the  parameters  of  motion  in  the  rest  frame  are 

,  3 

2Ttmy 

eB 


(2.24) 


l 


2  nn£  y“ 


eB‘ 


ixyd  . 
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Using  the  initial  conditions  of  Equation  (2.18),  we  obtain  the  analytical 
results, 


t  =  1 .7?x10"iZ  s 
d  =  8.53x10" 7  m 
l  =  2.68x10"S  m 


(2.25: 


Figure  2.4  presents  a  trajectory  plot  from  the  simulation.  Note 
that  the  cycle  appears  to  have  repeated  precisely  three  times  in  a 
simulation  of  duration,  3t.  The  measured  turning  distance  and  drift  are 


d  =  8 . 5x10" 7  m 
i  =  2.7x10" 6  m 


(2.26) 


These  values  were  obtained  by  graphical  measurement  of  the  trajectory 
plots. 


Table  2.3.  Input  data  for  Problem  2.3. 


title  *problem  2.3*  / 

comment  *drift  (e  cross  b)  in  constant  electromagnetic  field*  / 
system  1  / 

xlgrid  132  -4.81e-7  1  9.0e-6  9.0e-6  / 
x2grid  132  -4.07e-6  1  9.0e-6  9.0e-6  / 
conductor  boundary  +12232332322/ 
fields  1  1  300  1.78638e-14  / 
courant  00/ 
kinematics  10000/ 

forces  1  1  2  2  e2  null  null  6  b3  null  null  / 
function  e2  1  0  -3.0e7  / 
function  b3  1  0  -20.0  / 

populate  electrons  00. 012233  1.6e-19 

0.0534  0.4526  0.0  0.0  0.0  / 
diagnose  prscrib  0  2 

kinemat  0  301  / 

trajectory  300  300  1  -4.81e-7  8.52e-6  -4.07e-6  4.93e-6  / 
timeout  -10/ 
start  / 
stop  / 
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RGIC  VERSION:  SEPTEMBER  iS83  DATE :  85/01/1 
SIMULATION:  PROBLEM  2.3 

TRAJECTORY  PLOT  OF  ELECTRONS  t  ISPE  -  1  ) 

AT  TIME:  5.36E-12  SEC  FOR  300  TIME  STEPS 


Figure  2.4.  Trajectory  plot.  Problem  2. 
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DRIFT  (8  x  VB). 


2.4.1  Problem  Description. 

Spatial  variations  in  the  magnetic  field  will  result  in  a 
variety  of  particle  drifts.  One  drift  of  interest  in  plasma  physics  is 
that  due  to  a  transverse  gradient  in  the  magnetic  field.  Specifically,  we 
consider  the  fields, 

E  =  0 

B  =  zB  ,  (2.27) 

where  B  =  -  . 

y 

The  particle  initial  conditions  are 
t  =  x  =  z  -  0 

y  =  10" 2  m  (2.28) 

P  =  xPq  , 

where  p  =  9  x  10  m/sec. 

o 

2.4.2  Suggested  Approach. 

The  field  prescription  option  in  the  FORCES  command  along  with 
the  appropriate  FUNCTION  command  can  easily  be  used  to  provide  the  1/y 
dependence  in  the  magnetic  field  (see  function  no.  5  in  the  User's 
Manual).  Input  data  for  this  simulation  is  shown  in  Table  2.4. 
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2.4.3 


Analytical  Solution. 


This  problem  can  also  be  solved  analytically ,  at  least  for  the 
case  of  weak  spatial  gradients!  The  solution  involves  expansion  of  the 
magnetic  field  (or  gyromagnetic  frequency)  about  the  slowly  moving  orbital 
center,  or 


w(r)  =  -  — —  B(r)  =  a)  f  1  +  (—  —  )  y  •  r 

ymc  o  -  a  5y 

The  transverse  velocity  is  written  in  two  parts, 


(2.29) 


V  =  V  +  V  .  , 

o  d 


(2.30) 


where  v0  is  the  gyromagnetic  velocity  and  vp  is  a  small  drift 
velocity.  Substitution  into  the  frequency  expansion  yields  the  result, 


dv 


■  [  V  •  V  <y  •  r  >(  i-2  )]x»0. 
dt  B  by 


(2.31) 


Thus,  the  drift  velocity  has  a  nonvanishing  average  value, 


<vd>  =  ^  w0  *  <r0(y  *  ro)>  ’ 

B  oy 

or 


(2.32) 


<v  .>  =  14  (B  x  V  x  B)  . 
d  2B2 


(2.33) 


Using  the  parameters  of  Equations  (2.27)  and  (2.28),  we  obtain 
the  analytical  value  of  the  drift  velocity, 


t  3.  0.  Oackson,  Classical  Electrodynamics,  Bohn  Wiley  and  Sons,  New  York, 
1967  (p.  415). 
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Table  2.4.  Input  data  for  Problem  2.4. 


title  *problem  2.4*  / 

comment  *drift  (b  cross  delb)  in  spatially  varying 
magnetic  field*  / 
system  1  / 

xlgrid  132  -0.00026  1  0.0006  0.0006  / 

x2grid  132  0.00945  1  0.0006  0.0006  / 

fields  1  1  500  5.65233e-14  / 

conductor  boundary  +12232332322/ 

courant  00/ 

kinematics  10000/ 

forces  1116  null  null  b3x2  / 

function  b3x2  5  -1  -0.2  / 

populate  electrons  00. 012233  1.6e-19 

0.4339  0.9167  9 . 0e8  0.0  0.0  / 
diagnose  kinemat  0  501 

prscrib  0  501  / 

trajectory  500  500  1  -0.00026  0.00034  0.00945  0.01005  / 
timeout  -10/ 
start  / 
stop  / 
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FIELD  EMISSION. 


2.5.1  Problem  Description. 


Two  infinite  plates  separated  by  distance,  d,  are  maintained  at 
a  constant  potential  difference,  V.  We  arbitrarily  use  the  parameters, 


d  =  3  x  10-2  m 


(2.36) 


V  =  1  x  103  V 


The  cathode  is  assumed  to  break  down;  i.e.,  it  continuously  emits 
electrons  which  are  received  at  the  anode.  We  want  to  measure  the  final, 
steady-state  voltage  as  a  function  of  distance  and  the  incident  (cathode 
surface)  current  density  as  a  function  of  time. 


2.5.2  Suggested  Approach. 


The  one-dimensional  simulations  are  the  first  self-consistent 


problems  (i.e.,  mutual  interaction  between  particles  and  fields)  in  the 


set.  The  MAGIC  code  can  be  run  in  the  one-dimensional  mode  in  several 


ways.  We  suggest  using  a  single  cell  in  the  transverse  dimension  bounded 
by  periodic  symmetry  planes  (see  SYMMETRY  command).  The  periodic  symmetry 
model  requires  overlapping  cells  at  each  boundary;  thus,  four  full-grid 
points  (three  cells)  are  needed  in  the  transverse  dimension.  In  the 
longitudinal  direction,  any  number  of  cells  can  be  used;  however,  the 
quality  of  the  simulation  will  improve  with  careful  resolution  of  any 
space-charge  barrier  which  may  result. 


The  particles  for  this  simulation  must  be  introduced  using  a 
PARTICLES  command.  We  suggest  starting  the  particles  out  a  small  distance 
(fraction  of  the  first  cell)  from  the  cathode  and  with  a  commensurate 
velocity  (i.e.,  one  given  by  the  Child-Langmuir  equations).  These  are  all 
required  input  data.  For  output,  we  suggest  phase  space  plots  (p^  vs. 
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xi),  range  plots  of  fields  (E^  vs-  x  j.  and  3^  vs.  x^),  and  time  histories 
of  the  same  at  selected  locations. 

The  constant  voltage  for  this  simulation  can  be  applied  and 
maintained  using  the  CIRCUIT  command.  Note  that  the  abrupt  application  of 
a  constant  voltage  and  the  subsequent  breakdown  will  result  in  a 
transient.  That  is,  oscillations  will  occur  prior  to  reaching  a  steady 
state.  This  is,  of  course,  physical.  However,  it  has  been  shown^  that 
the  transient  can  be  minimized  if  the  voltage  is  initially  given  the 


time-dependence, 


where 


V(t)  =  V[  -  (  -) 
3  x 


.2  l/2 

T =  3(  OL.  ) 

2eV 


-  (-  )  ]  , 
3  x 


0  <  t  <  T 


(2.37) 


The  voltage  is  maintained  constant  at  V  after  the  period  t. 

The  simplest  way  to  input  this  function  is  to  use  the  numerical 
data  option.  The  input  data  for  this  simulation  is  shown  in  Table  2.5. 

2.5.3  Analytical  Solution. 


I 


In  one  dimension,  the  continuity  equation  is 


5  3  +  3.  p  =  0 
x  t 


(2.38) 


The  steady-state  assumption  requires  that 


T  M.  Lampel  and  M.  Tiefenback,  "An  Applied  Voltage  to  Eliminate  Current 
Transients  in  a  One-Dimensional  Diode,"  Appl.  Phys.  Lett.  43  (1),  3uly 
1983. 
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V  =  0 


and  subsequently, 


3  =  3+  , 


where  3+  is  the  current  density  at  the  cathode  surface, 
dimensional  Poisson  equation, 


=  -  - 
x  e 


(2.41  ) 


also  applies.  Here  the  density  is  given  by 


3 

P  =  V 


(2.42) 


The  velocity  can  be  related  to  potential  by  the  (nonrelativistic) 
expression, 

2e4>  i1/2 


v  =  (  21  ) 
m 


(2.43) 


Substitution  of  these  results  into  the  Poisson  equation  yields  the  result, 

.1/2 


a2  ♦  =  -  £  (  JH  )< 

e  2e 


(2.44) 


This  can  be  integrated  by  making  use  of  the  substitution, 

2 


a.  (  )  =  2(  «  H  £• )  , 

dx  dx  dx  dxz 


d  2  <t> 


(2.45) 
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to  obtain  the  familiar  result1", 


^  =  ( 


3 

2 


2/3 


) 

2e 


i/3 


.4/  3 


(2.46) 


Thus,  the  cathode  current  density  is  determined  by  anode  potential  and 
separation,  or 


9  m 

Using  the  parameters  in  Equation  (2.36),  we  obtain  the  analytical  result, 

3+  =  8.20x 104  A/m2  .  (2.48) 

Range  plot  results  from  the  simulation  are  illustrated  in 
Figures  2.6  and  2.7.  The  current  density  versus  distance  shown  in  Figure 
2.6  appears  to  be  stable  at  a  value  of 

3*  =  7.99x104  A/m2  .  (2.49) 

(The  slight  discrepancy  at  the  origin  is  due  to  the  way  local  charge 
conservation  is  handled  by  the  field  emission  algorithm.)  The  electric 
field  vs.  distance  is  shown  in  Figure  2.7;  this  result  exhibits  the 
expected  x1/3  behavior.  Finally,  Figure  2.8  presents  a  plot  of  particle 
phase  space.  For  this  nonrelativistlc  simulation,  the  particle  momentum 
should  go  as  x‘<  (i.e.,  the  square  root  of  potential). 


1/2  y  3/  2 


(2.47) 


t  G.  D.  Child,  Physical  Review  (Ser.  I)  32,  492  (1911). 
I.  Langmuir,  Physical  Review  2J_,  419  (  1923). 
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Table  2.5.  Input  data  for  Problem  2.5. 


title  *problem  2.5*  / 

comment  *space-charge  limiting  (child-langmuir  problem) *  / 
system  1  / 

xlgrid  1  42  2  0.0  40  1.875e-4  3 . 0e-2  / 
x2grid  1  4  2  0.0  2  3.75e-3  7.5e-3  / 
conductor  k  1  22  23/ 

conductor  a  -1  42  2  42  3  / 

symmetry  periodic  1  2  2  42  2  -1  2  3  42  3  / 
fields  1  1  4000  5.0e-13  / 
courant  00/ 

particles  emission  null  electrons  0. 00. 011221 

1 . 0e6  1 . 1680316e-5  0.0  0.0  0.0  0.0  0.0  12223/ 
forces  0.0  0.5  1.0  1.0  / 
poisson  2  a  +1.0  k  0.0  1  0.0  1.0  1  / 
circuit  voltage  10. 010  42  222/ 

function  voltage  0  10  0.0  0.0  6.0e-ll  1 . 67e4  1.2e-10  3.32e4 

1 . 8e-10  4 . 94e4  2.4e-10  6.46e4  3.0e-10  7 . 83e4 
3 . 6e-10  8 . 95e4  4.2e-10  9.71e4  4.8e-10  1 . 00e5 

1.0  1.00e5  / 

statistics  400  / 

diagnose  spacing  01/ 

linprint  4000  122  42  211/ 

linprint  4000  11  2  2  42  2  1  1  / 

trajectory  800  1  1  0.0  3.0e-2  0.0  3 . Oe-2  / 

phasespace  800  31  0  0.0  3.  Oe-2  0.0  3.75e-3  1.0e6  2 . 0e8 

-2.0e8  2 . 0e8  -2 . 0e8  2 . 0e8  / 
observe  1  0  twod  10.0  1.0  0.0  1.0  1  22  22 

twod  1  0.0  1.0  0.0  1.0  1  22  2  22  2 

twod  1  0.0  1.0  0.0  1.0  1  41  2  41  2  / 

range  800  1  722  42  23/ 

range  800  1  122  42  23/ 

range  800  1  11  2  2  42  2  3  / 

range  800  1  10  2  2  42  2  3  / 

timeout  -10/ 
start  / 


MAGIC  VERSION:  SEPTEMBER  1S83  DOTE:  34/12/13 
SIMULATION:  PROBLEM  2.5 


RANGE  PLOT  AT  TIME:  2.QOE-0S  SEC 
El  COMPONENT 

RANGING  FROM  (2,2)  TO  (42,2) 


C.OCu  G.3G5 


Figure  2.7. 
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2.6 


MAGNETIC  INSULATION. 


2.6.1 


Problem  Description. 


By  adding  a  magnetic  field  perpendicular  to  the  plane  described 
in  Problem  2.5,  a  one-dimensional  magnetic  insulation  simulation  is 
created.  For  example,  if  x  measures  distance  across  a  gap  of  width,  d, 
and  y  measures  distance  along  the  conducting  surfaces,  then  we  add  a 
constant  magnetic  field, 


B  =  zB 


(2.50) 


where  B  =  4.6  x  10"  T.  This  value  of  magnetic  field  was  selected  on  the 
basis  of  reasonable  hub  thickness.  Other  parameters  are  maintained  as  in 
Problem  2.5,  i.e. , 


d  =  3  x  10'2m 


V  =  1  x  10 3  V 


(2.51) 


As  with  the  Child-Langmuir  physics,  the  perpendicular  electric 
field,  Ex,  draws  electrons  out  of  the  cathode.  Now,  however,  the 
magnetic  field  tends  to  turn  electrons  in  the  y  direction.  This  leads  in 
turn  to  a  y-component  of  current  density,  a  y-component  of  electric  field 
and  a  z-component  of  (self)  magnetic  field.  Despite  the  presence  of  these 
nonvanishing  and  spatially  dependent  fields  (Ex,  Ev,  Bz,  3X, 
3y),  the  problem  is  still  one-dimensional. 

2.6.2  Suggested  Approach. 

The  external  magnetic  field  can  be  added  using  a  BEXTERNAL 
command.  (Note  that  the  simulation  automat ically  adds  this  field  to  the 
dynamic  field  from  Maxwell's  eguatlons  to  compute  particle  forces.) 
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The  time  transient  experienced  in  this  problem  is  even  more 
severe  than  that  in  Problem  2.5.  Therefore,  we  recommend  using  a  time- 
dependent  voltage  similar  to  that  suggested  in  2.5,  but  with  a  value  of  x 
equal  to  some  multiple  (say  4)  of  the  gyromagnetic  period.  Note  that  this 
period  corresponds  roughly  to  the  particle  lifetime  (surface-to-surface), 
as  shown  in  Problem  2.3.  Although  both  results  (stepped  or  transient 
voltage)  are  physical,  the  transient  case  corresponds  more  closely  to  what 
can  be  accomplished  with  actual  devices.  That  is,  a  slowly  increasing 
voltage  tends  to  mask  the  violence  of  the  boundary  layer  transient. 

The  input  data  for  this  simulation  is  shown  in  Table  2.6. 

2.6.3  Analytical  Solution. 

A  solution  of  this  problem  was  first  done  by  Sudan  and 
Lovelace*.  We  consider  that  the  problem  is  infinite  and  symmetrical  in 
the  y  and  z  directions,  so  that  all  quantities  will  be  unvarying  in  y  and 
z.  In  the  steady-state  solution  of  the  problem,  we  will  have  two 
regions,  one  near  the  cathode  filled  with  space  charge  and  one  next  to  the 
anode  which  is  a  vacuum  region.  The  regions  will  be  separated  by  a 
surface  at  x  =  xt. 

The  electron  equations  of  motion  are 


iv  =l£[E(x)+vB] 
t  x  m  L  x  y  z J 


(2.52) 


5  v  =  —  v  B 
t  y  m  y  z 


t  R.  N.  Sudan  and  R.  V.  Lovelace,  Phys.  Rev.  Lett.  31,  1174,  (1973). 


2-32 


urw  v  v  nn  'J<  vv/v .v .'v»  .-juvAr  v  .vw.th  tt  ■/■  -- 


Since  we  assume  that  vx  =  vy  =  0  for  electrons  emerging  from  the 

cathode  (x  =  0)  and  that  Bz  is  constant  in  space,  we  can  integrate  both 
sides  of  Equation  (2.52)  to  obtain 


V  =  XCJ 

y  c 


(2.53) 


where  wc  =  eB0/m  is  the  electron  cyclotron  frequency.  We  set  M  0  at 
the  cathode  and  obtain  through  energy  conservation  the  equation 


—  (v2  +  v2)  =  e  <t>  ( x ) 

■)  x  y 


(2.54) 


This  requires  that,  at  the  turning  point, 


1  2  2  e  . 

-  “  x  =  _  4»(x  ) 

2  c  t  m  t 


(2.55) 


We  assume  that  there  will  be  a  flow  of  charge  from  the  cathode  in  this 
problem  that  is  balanced  by  an  opposing  return  flow.  Both  flows  will 
be  separately  conserved,  except  at  the  turning  point.  We  therefore  have 
by  conservation  of  charge: 


5  3  +  53  =5.  p  , 

xx  y  y  t 


(2.56) 


where  3X  and  3y  and  the  current  density  in  the  x  and  y  direction 
respectively  and  p  is  the  change  density.  Since  we  have  steady-state 
conditions  5tp  =  0  and  also  5y  3y  =  0,  because  all  quantities  must 

be  constant  in  the  y  direction,  3X  =  30  =  pvx  =  constant.  Because 

of  both  emitted  and  return  flows,  the  density  must  be 


(2.57) 
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We  have  both  magnetic  and  electric  fields  in  the  problem.  We  approximate 
the  magnetic  fields  as  being  constant  in  space,  Bz  =  B0.  The  electric 
fields  are  found  via  the  electrostatic  potential  as  determined  from 
Poisson's  equation: 


d2<t>  =  -  —  . 
x  e 

o 


(2.58) 


For  the  vacuum  region,  p  =  0,  so  <t>  will  be 


$  =  V  +  f(d-x)  , 
o 


(2.59) 


where  f  is  a  constant  to  be  determined  as  part  of  the  solution  of  the 
problem.  The  potentials  and  derivatives  for  the  two  regions  will  be 
matched  at  xt. 

The  cathode  surface,  at  steady- state,  is  a  space  charge  limited 
emitter,  implying 


9  4>(0)  =  0  . 


By  Poisson's  equation  and  conservation  of  charge  we  have 


a**  =  -  2 ]/ev. 

X  O  O  X 

We  can  use  energy  conservation  and  the  equations  of  motion  to  write 


d2<t  =  (m/2e)  1/2  3  /(  *  -  3_  x2)  l/  2 
x  0  2 


.  2  m  2 

where  q  =  —  w 
e  c 


(2.60) 


(2.61 ) 


(2.62) 
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We  define 


2  2 

4i  =  0  -  _  qz  xz 


(2.63: 


and  write  (using  prime  notation,  5X tp  =  4-' ) 


.„  2  rK  >1,2 

^  +  q  =  C-jp—  J  / 


where  K  =  23  (m/2e) 
o 


This  can  be  integrated  to  obtain 
4»'  =  (4K  4-1/2  -  2q2  4>)l/2  . 


(2.64) 


(2.65) 


This  equation  gives  us  a  condition  on  4>'  and  thus  0'  at  the  turning 
point.  The  variable  4>  is  actually  (m/2e)vx2.  Therefore,  0  =  0  at  xt 


so  that 


'(xt)  =  q  xt  . 


(2.66) 


Using  the  fact  that  0  and  0'  must  be  continuous  across  the  surface  at 
xt,  we  have  then  the  condition  for  the  turning  point, 


=  -L-  B2x  (d  -  x,/2)  . 

o  2  o  t  t 

me 


(2.67) 


The  steady-state  field  in  the  gap  is  then  also  the  same  as  for 
parapotent ial  flow, 


(2.68) 
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For  the  parameters  in  Equation  (2.51),  we  have 


vV. 

v'-JV 


=  1.18x10“  m 


E  =  4. 14x10°  M/m 


(2.69) 


Results  from  the  simulation  are  shown  in  Figures  2.9  through 
2.12.  The  transverse  current  density  2y  is  shown  in  Figure  2.9.  The 
normal  electric  field  Ex  as  a  function  of  distance  is  shown  in  Figure 
2.10.  From  this  figure,  we  estimate  the  turning  distance  and  gap  field  to 


x^  =  1.2x10“  m 


E  =  4.2x10  V/m  , 


(2.70) 


in  reasonable  agreement  with  the  analytical  results.  Figure  2.11 
illustrates  the  seif-magnetic  field  (dynamic  component  only )  vs. 
distance.  The  diamagnetic  effect  is  confirmed  small  in  comparison  to  the 
externally  applied  field,  thus  justifying  the  initial  assumption. 
Finally,  Figure  2.12  presents  a  particle  phase-space  ( p 2  vs.  xj)  plot. 
(This  plot  is  composed  of  individual  particle  coordinates.)  Note  the 
extent  to  which  this  result  agrees  with  Equation  (2.53).  Graphically 
taking  the  slope  of  this  curve,  we  estimate  from  the  simulation  to  be 
9x10  rad/sec  vs.  the  analytical  value  of  9.07x10  rad/sec. 
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Table  2.6.  Input  data  for  Problem  2.6. 


title  *problem  2 . 6i*  / 

comment  ^magnetic  insulation  boundary  layer*  / 
system  1  / 

xlgrid  1  42  2  0.0  40  1.875e-4  3.0e-2  / 
x2grid  1  4  2  0.0  2  3.75e-3  7.5e-3  / 
conductor  k  1  22  23/ 

conductor  a  -1  42  2  42  3  / 

symmetry  periodic  1  2  2  42  2  -1  2  3  42  3  / 
fields  1  1  12C00  5.0e-13  / 
courant  00/ 

particles  emission  null  electrons  0.00.011221 

1 . 0e6  1 . 1680316e-5  0.0  0.0  0.0  0.0  0.0  1  2  2  2  3  / 
forces  0.5  1.0  1.0  / 
bexternal  0.0  0.0  4.6176342e-2  / 
poisson  2  a  +1.0  k  0.0  1  0.0  1.0  1  / 
circuit  voltage  10. 010  42  222/ 

function  voltage  0  17  0.0  0.0  2.065e-10  8.896e3  4.13e-10  1.778e4 
6 . 195e-10  2 . 664e4  8.26e-10  3.542e4  1.0325e-09  4.41e4 
1.239e-09  5.25e4  1.4455e-09  6.069e4  1.652e-09  6.85e4 
1 . 8585e-09  7.57e4  2.065e-09  8.24e4  2.27l5e-09  8 . 82e4 
2 . 478e-09  9 . 31e4  2.6845e-09  9 . 68e4  2.891e-09  9 . 92e4 

3.0975e-09  1.00e5  1.0  1.00e5  / 

statistics  400  / 
diagnose  601 

3  4000  4001 

4  4000  4001 

5  4000  4001  / 

linprint  4000  122  42  211/ 

linprint  4000  11  2  2  42  2  1  1  / 

trajectory  800  1  1  0.0  3.0e-2  0.0  3.0e-2  / 
phasespace  800  31  0  0.0  3.0e-2  -1.0  1.0  -1.0e8  1 . 0e8 

-9 . 0e9  9 . 0e9  -9.0e9  9.0e9  / 
phasespace  800  41  0  0.0  3.0e-2  -1.0  1.0  -9.0e9  9 . 0e9 

-1.0e8  1 . 0e8  -9.0e9  9.0e9  / 
observe  1  twod  10.01.01  22  22 

twod  10. 01. 01  22  2  22  2 
twod  10.01.01412412 
twod  10. 01. 06  22  22 

twod  10. 01. 06  22  2  22  2 


i 

E 

Table  2.6 

• 

Tnput  data  for 

Problem  2.6.  (continued) 

$ 

m 

\ 

twod 

1 

0.0  1.0  6 

41 

2  41 

2 

f, 

twod 

1 

0.0  1.0  8 

3 

2  3 

2 

i 

twod 

1 

0.0  1.0  8 

22 

2  22 

2 

I 

twod 

1 

0.0  1.0  8 

41 

2  41 

2 

"J 

twod 

1 

0.0  1.0  7 

2 

2  2 

2 

“j 

twod 

1 

0.0  1.0  7 

22 

2  22 

2 

it 

twod 

1 

0.0  1.0  7 

41 

2  41 

2 

c 

twod 

1 

0.0  1.0  10 

2 

2  2 

2 

1 

twod 

1 

0.0  1.0  10 

22  2  22 

2 

twod 

1 

0.0  1.0  10 

41  2  41 

2  / 

j  range 

800 

1 

7  2 

2 

42  2  3  / 

l*1  .  range 

800 

1 

8  2 

2 

42  2  3  / 

*4  range 

800 

1 

8  2 

3 

42  3  3  / 

range 

800 

1 

1  2 

2 

42  2  3  / 

[,  range 

800 

1 

1  2 

3 

42  3  3  / 

-•  range 

800 

1 

2  2 

2 

42  2  3  / 

I-  range 

800 

1 

2  2 

3 

42  3  3  / 

l’  range 

800 

1 

6  2 

2 

41  2  3  / 

j  range 

800 

1 

6  2 

3 

41  3  3  / 

range 

800 

1 

11  2 

2 

42  2  3  / 

/•,  timeout  60 

1  / 

V  start 

/ 

SECTION  3 

PULSED  POWER  TRANSPORT 


This  section  presents  problems  involving  pulsed  power  transport 
(e.g.,  diode)  applications.  The  objective  is  to  combine  previously  intro¬ 
duced  simulation  concepts  in  realistic  applications.  However,  this  objec¬ 
tive  must  be  tempered  by  economic  reality.  Realistic  applications  will 
often  require  expensive  simulations,  whereas  the  simulations  presented 
here  must  be  sufficiently  inexpensive  to  allow  use  as  a  learning  tool. 
Thus,  we  have  established  a  maximum  cost  of  about  five  minutes  of  Cray  CPU 
time  per  problem,  and  have  compromised  simulation  parameters  to  meet  this 
restriction.  The  resulting  problems  may  lack  geometric  and  temporal 
realism;  however,  they  introduce  the  important  simulation  concepts  and 
are,  at  the  same  time,  affordable. 


i 


3.1  COAXIAL  LINE  WITH  FIELD  EMISSION. 

3.1.1  Problem  Description. 

The  first  problem  involves  a  uniform  section  of  coaxial  trans¬ 
mission  line,  as  previously  illustrated  in  Figure  1.1.  For  convenience, 
the  dimensions  of  Problem  1.1, 

r0  =  0.20m 

rj  =  0.10m  (3.1) 

L  =  0.60m  , 

are  duplicated  exactly.  Whereas  Problem  1.1  simulated  a  cold-test,  here 
we  wish  to  cause  field  emission  of  electrons  to  occur  on  the  cathode. 
Therefore,  the  incident  voltage  pulse  Is  semi-infinite  with  a  4  MV  peak 
and  a  rise  time  of  1  nsec,  or 

wi)  ,  o  <  t  <  t 

p  T  ; 

0(t)  =  (3.2) 

4>+  ,  x  <  t  <  ®  . 

P 

Plasma  formation  will  be  allowed  on  the  cathode  at  a  breakdown  field  of 
2.5x107  V/m. 

We  wish  to  measure  the  voltage  and  magnetic  field  (current)  at 
the  inlet,  midpoint,  and  outlet  as  a  function  of  time,  the  radial  electric 
field  and  axial  current  density  as  a  function  of  radius,  and  the  particle 
trajectories  at  selected  times  during  the  simulation.  The  primary 
objective  is  to  observe  the  following:  (1)  field  emission  at  the  cathode 
and  electron  leakage  to  the  anode  as  the  incident  wave  front  propagates 
down  the  line,  (2)  magnetic  self-insulation  behind  the  incident  wave 
front,  and  (3)  the  effect  of  the  space-charge  boundary  layer  on  the 
impedance  of  the  line. 
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3.1.2  Suggested  Approach. 

Resolution  of  the  space-charge  barrier  in  field  emission 
typically  requires  modification  to  the  spatial  grid.  In  Problem  1.1,  it 
was  possible  to  use  a  uniform  spatial  grid  in  both  dimensions.  However, 
the  particle  simulations  of  Problems  2.5  and  2.6  show  that  non-uniform 
spacing  normal  to  an  emitting  surface  produces  more  accurate  results 
due  to  better  resolution  of  gradients.  In  MAGIC,  non-uniform  regions  can 
be  defined  over  portions  of  the  grid  and  matched  to  other  non-uniform 
regions  or  to  uniform  regions.  It  is  important  to  make  sure  that  when 
matching  regions,  the  cells  bordering  an  adjacent  region  should  have 
approximately  the  same  size.  Within  a  region,  a  20%  cell-to-cell 
variation  is  generally  acceptable.  Also,  the  cell  aspect  ratio  should  not 
exceed  five.  In  this  case,  we  have  used  a  spacing  of  0.002  m  at  the 
cathode,  increasing  linearly  to  0.008  m  at  the  anode. 

All  boundary  conditions  are  the  same  as  in  Problem  1.4. 
However,  the  implicit  (time-biased)  field  algorithm  is  recommended  to 
stabilize  the  representation  of  field  emission.  Two-dimensional  particle 
kinematics  is  sufficient  for  this  simulation  (and  all  others  in  this  set), 
since  the  initial  distribution  is  two-dimensional  and  the  TM  mode  cannot 
induce  a  three-dimensional  component. 


Finally,  electron  field  emission  must  be  specified  along  the 
cathode  surface.  In  MAGIC,  field  emission  must  be  enabled  for  each 
intended  (conformal)  surface  and  species  using  the  PARTICLES  command. 
Note  that  field  emission  can  occur  (depending  on  the  dynamics)  only  on 
surfaces  which  have  been  so  enabled.  Also,  each  command  specifies  only 
one  particle  species  (e.g.,  electrons).  Thus,  if  emission  can  occur  at 
one  surface  for  more  than  one  species,  each  must  be  enabled  separately. 
(A  common  example  Involves  electron  beam  emission  and  ion  field  emission 
from  the  same  surface.)  In  this  case,  we  also  recommend  leaving  a  five 
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cell  zone  at  the  inlet  in  which  no  emission  occurs  to  allow  the  TEM 
waveform  to  develop. 

This  input  data  for  this  simulation  is  shown  in  Table  3.1. 

3.1.3  Analytical  Solution. 


The  space-charge  boundary  layer  can  be  considered  as  a  region  of 
impedance  (Z2)  different  from  that  of  the  inlet  region  (Z^).  Then  the 
steady-state  relations  developed  previously  in  Equation  (1.15)  for  an 
impedance  mismatch  are  valid.  The  voltage  equation  and  the  measured 
simulation  voltage  can  then  be  used  to  estimate  the  boundary  layer  region 
impedance, 


12 


< t 


(3.3) 


From  the  inlet  voltage  shown  in  Figure  3.1,  we  estimate  the  final  voltage 
to  be  <(>  -  3.25  MV.  Then  Equation  (3.3)  yields 

Z2  -  28.5  Q  .  (3.4) 

Figures  3.2  and  3.3  illustrate  the  progression  of  the  pulse  down  the  line. 

The  same  process  can  be  repeated  using  the  current  equation.  In 
this  case,  the  magnetic  field  In  the  gap  (at  the  anode)  should  be  used  to 
compute  the  total  current,  or 

Io  =  ■TT0e<r°)  •  (3-5> 


If  the  magnetic  field  at  the  cathode  is  used,  then  the  current 
contribution  due  to  the  boundary  layer  must  be  included,  yielding 


1l  =  B0(ri)  + '  dr  2ur  Vr)  • 

r . 

1 


(3.6) 


In  principle,  the  two  currents  (Ij  and  Iq)  should  be  identical,  and 
either  can  be  used  in  Equation  (1.15)  to  estimate  the  impedance. 

The  impedance  can  also  be  estimated  more  crudely  by  measuring 
the  hub  thickness  (or  effective  gap)  and  using  the  conventional  impedance 
equations.  From  the  particle  trajectory  plot  shown  in  Figure  3.4,  we 
estimate  the  effective  cathode  radius  (near  outlet)  to  be  r  =  0.15  m. 
Then  Equation  (1.9)  yields 


Z2  =  17.3  Q 


(3.7) 


The  problem  with  this  commonly  used  method  is  that  it  ignores  the  boundary 
layer  dynamics. 


Table  3.1.  Input  data  for  Problem  3.1. 


title  ^problem  3.1*  / 

comment  *base  problem:  41.6  ohm  coaxial  transmission  line.*  / 
system  2  / 

xlgrid  1  62  2  0.0  60  0.01  0.6  / 
x2grid  1  22  2  0.1  20  0.002  0.1  / 
fields  1  3  800  6.25e-12  0.5  0.5  0.0 

4  1.0  0.29912  0.15022  0.11111  / 
courant  00/ 

conductor  cathode  12  2  62  2  / 
conductor  anode  -1  2  22  62  22  / 
particles  emission  null  electrons  0061121 

1.0e+7  6 . 25e-5  2.5e7  0.0  0.0  0.0  0.0  1  7  2  62  2  / 

kinematics  10000/ 
currents  64  1.0  / 
forces  0  0.5  1.0  1.0  / 

voltage  incident  radial  1  twod  01. 01222  22  2/ 

function  incident  0  3  0.0  0.0  1.0e-9  4.0e6  1.0  4.0e6  / 

function  radial  5  -11/ 

lookback  1  1  -1  62  2  62  22  / 

diagnose  spacing  0  1  courant  01/ 

statistics  200  / 

trajectory  200  1  1  0.0  0.6  -0.15  0.45  / 
observe  1  1  twod  10. 01. 00. 01. 02  2  2  2  22 

twod  1  0. 0  1. 0  0. 0  1. 0  2  32  2  32  22 

twod  10. 01. 00. 01. 02  62  2  62  22 

twod  10.01.00.01.06  2  2  2  2 

twod  10. 01. 00. 01. 06  32  2  32  2 

twod  10. 01. 00. 01. 06  61  2  61  2 

twod  1  0.0  1.0  0.0  1.0  6  2  21  2  21 

twod  1  0.0  1.0  0.0  1.0  6  32  21  32  21 

twod  1  0.0  1.0  0.0  1.0  6  61  21  61  21  / 

range  200  12  22  2  22  1 

1  2  32  2  32  22  1 

1  2  62  2  62  22  1 

16  22  2  22  1 
1  6  32  2  32  22  1 

1  6  61  2  61  22  1  / 

display  0  0.0  0.6  -0.15  0.45  / 
output  0  / 
timeout  60  1  / 
start  / 
stop  / 
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SIMULATION:  PROBLEM  3.1 
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Figure  3.3.  Outlet  voltage  vs.  time,  Problem  3.1 
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COAXIAL  LINE  WITH  SHORT. 


3.2.1  Problem  Description. 

Problem  3.2  uses  the  coaxial  transmission  line  of  Problem  3.1, 

but  with  a  short  across  the  outlet.  Repeat  the  measurements  defined 
above . 

3.2.2  Suggested  Approach. 

Replace  the  LOOKBACK  command  at  the  outlet  with  a  CONDUCTOR 
command,  which  will  short  the  cathode  and  anode.  Input  data  for  this 

simulation  is  shown  in  Table  3.2. 

3.2.3  Analytical  Solution. 

Problem  3.2  is  similar  to  the  classical  three-region  impedance 
problem,  illustrated  in  Figure  3.5.  In  this  case,  Z^  represents  the  inlet 
(plasma-free)  region,  Z2  represents  the  boundary  layer  (plasma)  region, 
and  Z3(=0)  represents  the  short.  The  present  analysis  is  complicated  by 
the  fact  that  the  boundary  layer  effect  is  transient;  that  is,  the  value 
of  Z2  changes  during  the  period  in  which  it  has  an  effect.  This  is  in 
contrast  with  Problem  3.1,  in  which  the  boundary  layer  achieves  a 

non-trlvial  ( Z 2  *  Z 1 )  steady  state.  The  following  derivation  assumes  a 
constant  value  of  Z 2  as  well  as  a  steady-state  pulse. 

From  our  previous  derivation  for  a  single  mismatch,  we  obtained 
the  scattered  wave  solution  ($*)  for  an  incident  wave  ($+)  from  the 

left-hand  side,  or 

<tr  =  ao*  ,  (3.8) 
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where  the  scattering  coefficient  is 


1 2  -  Zi 

Ti  +  z2 


(3.9) 


If  the  wave  is  incident  from  the  right-hand  side,  then  the  solution  is 


=  -a<t+  . 


(3.10) 


Figure  3.3  exhibits  two  mismatches  ( Z i , 2 2  and  Z2,Z3),  with 
scattering  coefficients  labeled  a  and  b,  respectively.  Thus,  multiple 
scatterings  occur,  as  waves  are  trapped  in  region  2  between  the 
mismatches.  The  solution  is  given  by  the  infinite  series, 

<t>/4>+  =  (1  +•  a)  +  (1  +  a ) b(  1  -  a) 

+  (1  ♦  a)b(1  -  a)(-ab)  +  ...  ,  (3.11) 

where  the  terms  represent  contributions  at  increasing  times,  n6t,  where  6t 
is  the  transit  time  in  region  2  (between  a  and  b).  The  series  may  be 
summed  to  yield 

0  =  <t>+  — - -Z-3—  ,  t—  ,  (3.12) 

Zi  +  Z3 

which  result  should  be  compared  with  Equation  M.15).  In  other  words,  the 
steady-state  response  is  independent  of  the  value  of  Z2.  In  some  cases, 
the  transient  response  may  be  highly  dependent  upon  Z2,  since  the  terms  in 
Equation  (3.11)  add  sequentially  during  the  transient. 

For  this  simulation,  the  transient  response  is  illustrated  in 

the  voltage  plots  of  Figures  3.6  and  3.7.  The  inlet  voltage  in  Figure  3.6 

closely  resembles  the  previous  result  in  Figure  3.1.  (The  scattered  wave 
from  the  short  is  just  beginning  to  reach  the  inlet.)  The  effect  of  the 

short  is  clear  in  Figure  3.7,  which  shows  the  voltage  at  midpoint.  (This 

result  should  be  comoared  with  that  of  Figure  3.2.) 
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One  important  point  is  that  impedance  effects  due  to  plasma  can 
be  strongly  influenced  by  effects  in  other  regions.  This  is  illustrated 
in  the  electron  trajectory  plots  in  Figures  3.S,  3.9,  and  3.10,  which 

represent  different  points  of  time  (2.50,  3.75,  and  5.00  nsec)  in  the 

simulation.  During  the  transient,  the  impedance  in  this  region  (Z2)  will 
change  from  41  .6  ohms  to  a  very  low  value  at  the  leading  edge  and  then 
from  approximately  29.5  ohms  to  a  final  value  of  41.6  ohms. 

It  is  also  useful  to  compare  the  responses  of  Problems  3.2  and 
1.3  (from  the  earlier  set).  This  illustrates  the  effect  of  the  plasma 
upon  propagation  velocity. 


.* 


Table  3.2.  Input  data  for  Problem  3.2. 


title  *problem  3.2*  / 

comment  *problem3.2  duplicates  3.1, 

but  with  a  short  at  the  outlet.*  / 
system  2  / 

xlgrid  1  62  2  0.0  60  0.01  0.6  / 
x2grid  1  22  2  0.1  20  0.002  0.1  / 
fields  1  3  800  6.25e-12  0.5  0.5  0.0 

4  1.0  0.29912  0.15022  0.11111  / 
courant  00/ 

conductor  cathode  12  2  62  2  / 

conductor  anode  -1  ?  22  62  22  / 
conductor  short  -1  62  2  62  22  / 
particles  emission  null  electrons  0061121 

1.0e+7  6 . 25e-5  2.5e7  0.0  0.0  0.0  0.0  1  7  2  62  2  / 
kinematics  10000/ 
currents  64  1.0  / 
forces  0  0.5  1.0  1.0  / 

voltage  incident  radial  1  twod  01.01222222/ 
function  incident  0  3  0.0  0.0  1.0e-9  4.0e6  1.0  4.0e6  / 
function  radial  5  -11/ 

z  lookback  1  1  -1  62  2  62  22  / 
diagnose  spacing  0  1  courant  01/ 
statistics  200  / 

trajectory  200  1  1  0.0  0.6  -0.15  0.45  / 
observe  1  1  twod  10. 01. 00. 01. 02  2  2  2  22 

twod  10. 01. 00. 01. 02  32  2  32  22 

twod  10. 01. 00. 01. 02  62  2  62  22 

twod  10.01.00.01.06  2  2  2  2 

twod  10. 01. 00. 01. 06  32  2  32  2 

twod  10. 01. 00. 01. 06  61  2  61  2 

twod  1  0.0  1.0  0.0  1.0  6  2  21  2  21 

twod  1  0.0  1.0  0.0  1.0  6  32  21  32  21 

twod  1  0.0  1.0  0.0  1.0  6  61  21  61  21  / 

range  200  12  22  2  22  1 

1  2  32  2  32  22  1 

1  2  62  2  62  22  1 

16  22  2  22  1 
1  6  32  2  32  22  1 

1  6  61  2  61  22  1  / 

display  0  0.0  0.6  -0.15  0.45  / 
output  0  / 
timeout  60  1  / 
start  / 
step  / 
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C  VERSION:  SEPTEMBER  1933  CATE : 
SIMULATION:  PROBLEM  3.2 

TRAJECTORY  PLOT  OF  ELECTRONS  (  ISPE 
AT  TIME:  2.50E-09  SEC  FOR  1  TIME  S 


re  3.8.  Trajectory  plot 
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COAXIAL  LINE  WITH  GEOMETRIC  MISMATCH. 


3.3.1  Problem  Description. 

The  previous  problem  illustrated  the  effect  of  a  short  (Z  =  0) 
on  the  boundary  layer  dynamics.  Here,  we  consider  a  geometric 

modification  of  Problem  3.1  in  which  the  cathode  radius  ( r ^ )  is  reduced 
to  0.05  m  at  a  distance  0.2  m  from  the  outlet.  In  effect,  this  doubles 
the  (cold-test)  impedance  of  the  outlet  region. 

3.3.2  Suggested  Approach. 

Non-uniform  spacing  is  now  suggested  for  the  axial  grid  as  well 
as  the  radial.  If  the  spacing  is  made  symmetric  about  the  mismatch,  then 
a  cell  aspect  of  unity  is  attained  at  the  corner  of  the  shoulder.  We 

suggest  limiting  the  field  emission  to  the  low  impedance  surface  (inlet 

region)  of  the  line.  This  input  data  for  this  simulation  is  shown  in 

Table  3.3. 

3.3.3  Analytical  Solution. 

Results  developed  for  Problem  3.2  are  applicable.  The  inlet 
voltage  is  shown  in  Figure  3.11.  From  our  measurement  of  the  voltage  (  t>  - 
3.25  MV)  and  Equation  (3.12),  we  estimate  the  outlet  impedance  to  be 

Z3  =  23.4  2  ,  (3.13) 

almost  the  same  as  that  for  Z2  in  Problem  3.1.  In  other  words,  the  outlet 
region  impedance  Is  dominated  by  field  emission  from  the  inlet  region. 

Figure  3.12  presents  an  electron  trajectory  plot  at  the  end  of 
the  simulation. 


Table  3.3.  Input  data  for  Problem  3.3 


title  ^problem  3.3*  / 

comment  *in  problem  3.3  an  impedance  mismatch  is  introduced 
20  cells  from  the  outlet.*  / 
system  2  / 

xlgrid  1  62  2  0.0  60  0.01  0.6  / 

x2grid  1  32  2  0.05  10  0.008  0.05  20  0.002  0.1  / 
fields  1  3  800  6.25e-12  0.5  0.5  0.0 

4  1.0  0.29912  0.15022  0.11111  / 
courant  00/ 

conductor  cathode  1  2  12  42  12  42  2  62  2  / 

conductor  anode  -1  2  32  62  32  / 

particles  emission  null  electrons  0061121 

1.0e+7  6 . 25e-5  2.5e7  0.0  0.0  0.0  0.0  1  7  12  42  12  / 
kinematics  10000/ 
currents  64  1.0  / 
forces  0  0.5  1.0  1.0  / 

voltage  incident  radial  1  twod  01. 012  12  2  32  2/ 

function  incident  0  3  0.0  0.0  1.0e-9  4.0e6  1.0  4.0e6  / 

function  radial  5  -11/ 

lookback  1  1  -1  62  2  62  32  / 

diagnose  spacing  0  1  courant  01/ 

statistics  200  / 

trajectory  200  1  1  0.0  0.6  -0.15  0.45  / 
observe  1  1  twod  1  0.0  1.0  0.0  1.0  2  2  12  2  32 

twod  1  0.0  1.0  0.0  1.0  2  32  12  32  32 

twod  1  0.0  1.0  0.0  1.0  2  42  12  42  32 

twod  10. 01. 00. 01. 02  62  2  62  32 

twod  1  0.0  1.0  0.0  1.0  6  2  12  2  12 

twod  1  0.0  1.0  0.0  1.0  6  32  12  32  12 

twod  1  0.0  1.0  0.0  1.0  6  42  12  42  12 

twod  10. 01. 00. 01. 06  61  2  61  2 

twod  1  0.0  1.0  0.0  1.0  6  2  31  2  31 

twod  1  0.0  1.0  0.0  1.0  6  32  31  32  31 

twod  1  0.0  1.0  0.0  1.0  6  42  31  42  31 

twod  1  0.0  1.0  0.0  1.0  6  61  31  61  31  / 


52  2  62  32 

2  12  2  12 


51  2 
2  31 


61  2 
2  31 


range  200  12  2  12  2  32  1 

1  2  32  12  32  32  1 

1  2  42  12  42  32  1 

1  2  62  2  62  32  1 

16  2  12  2  32  1 
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3.4 


COLD-TEST  OF  DIODE. 


3.4.1  Problem  Description. 

Based  on  the  coaxial  line  dimensions  in  Problem  3.3,  we  have 

created  a  diode  by  reducing  both  anode  and  cathode  radii  in  the  outlet 
region  until  they  extend  to  the  axis.  The  geometry  and  dimensions  are 
shown  in  Figure  3.13.  Mo  emission  surfaces  are  specified  in  this  problem, 
so  cold-test  results  are  expected. 

3.4.2  Suggested  Approach. 

The  spatial  grid  can  be  decreased  by  0.2  m  in  the  axial 
coordinate,  but  it  needs  to  be  increased  by  0.1  m  in  the  radial 
coordinate.  Because  the  problem  involves  cold-test  (no  particles),  a 
relatively  coarse  spacing  could  be  used.  However,  three  subsequent 
simulations  (Problems  3.5,  3.6,  and  3.9)  will  require  particles  for  the 
same  geometry.  Therefore,  non-uniform  spacing  similar  to  that  used  in 

Problem  3.3  is  recommended  so  that  the  electromagnetic  model  will  be 

identical  for  all  diode  simulations. 

Since  this  simulation  extends  down  to  zero  radius,  it  will  be 
necessary  to  use  the  SYMMETRY  command  to  specify  an  axial  symmetry 
boundary  condition.  (All  symmetry  boundary  conditions  in  MAGIC  (axial, 
mirror,  and  periodic)  must  be  specified  explicitly  by  the  user.  They  are 
not  supplied  automat ically  by  the  code.)  Mote  that  the  boundary  need  not 
extend  as  far  as  is  shown  in  Figure  3.13,  but  it  must  cover  the 
non-conducting  section  on  axis.  Time  histories  and  range  plots  should  be 
placed  at  the  inlet,  the  first  impedance  mismatch,  and  the  axis  of 

symmetry. 


Input  data  for  this  simulation  is  shown  in  Table  3.4. 


3.4.3 


Analytical  Solution. 


The  analytical  value  for  the  impedance  in  the  cylindrical  region 
between  cathode  and  anode  is 


Z 


1—  (  ^  l1/ 2  d 

2ti  v  £q  r 


(3.14) 


where  d  represents  the  gap  between  conductors  and  r  is  the  radius.  Once 
again,  the  analytical  solution  discussed  in  Problem  3.2  is  applicable.  In 
this  case,  the  impedance  I2  is  variable  according  to  Equation  (3.14),  and 
we  can  envision  a  third  region  of  infinite  impedance  (Z3  =  ®)  at  the 
axis.  Thus,  the  series  solution  developed  in  Problem  3.2  is  applicable 
for  the  transient  response.  This  will  be  most  apparent  in  the  measurement 
of  inlet  voltage  versus  time  shown  in  Figure  3.14,  where  several  of  the 
multiple  scatterings  described  in  Equation  (3.11)  can  be  observed. 
Finally,  according  to  Equation  (3.12),  the  final  steady-state  voltage 
should  be  exactly  twice  the  incident  voltage,  as  confirmed  in  Figures  3.15 
and  3.16. 
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Table  3.4.  Input  data  for  Problem  3.4. 


title  ^problem  3. 4*  / 

comment  ’‘problem  3.4  is  a  cold  test  of  a  diode.*  / 
system  2  / 

xlgrid  1  67  2  0.0  25  0.008  0.2  20  0.008  0.1 

20  0.002  0.1  / 

x2grid  1  42  2  0.0  20  0.008  0.1  20  0.002  0.1  / 
fields  1  3  800  6.25e-12  0.5  0.5  0.0 

4  1.0  0.29912  0.15022  0. 11111  / 
courant  00/ 

conductor  cathode  1  2  22  47  22  47  2  / 
conductor  anode  -1  2  42  67  42  67  2/ 
symmetry  axial  12  2  67  2  / 

z  particles  emission  null  electrons  0061121 

1.0e+7  6 . 25e-5  2.5e7  0.0  0.0  0.0  0.0  1  7  22  47  22  / 
kinematics  10000/ 
currents  64  1.0  / 
forces  0  0.5  1.0  1.0  / 

voltage  incident  radial  1  twod  01. 012  22  2  42  2/ 
function  incident  0  3  0.0  0.0  1.0e-9  4.0e6  1.0  4 . 0e6  / 
function  radial  5  -11/ 

diagnose  spacing  0  1  courant  01/ 
statistics  200  / 

trajectory  200  1  1  0.0  0.6  -0.15  0.45  / 
observe  1  1  twod  1  0.0  1.0  0.0  1.0  2  2  22  2  42 

twod  1  0.0  1.0  0.0  1.0  2  47  22  47  42 

twod  10. 01. 00. 01. 01  47  2  67  2/ 

range  200  12  2  22  2  42  1 

1  2  47  22  47  42  1 

1  1  47  2  67  21/ 

display  0  0.0  0.6  -0.15  0.45  / 
output  0  / 
timeout  60  1  / 
start  / 
stop  / 
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Figure  3.14.  Inlet  voltage  vs.  time,  Problem  3.4. 
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3.5.1 


OIODE  WITH  CATHODE  FIELD  EMISSION. 


Problem  Description. 


Problem  3.5  is  the  same  as  Problem  3. A,  except  that  electron 
field  emission  is  enabled  everywhere  on  the  cathode  (two  separate,  confor¬ 
mal  surfaces).  The  field  emission  and  incident  pulse  parameters  are  the 
same  as  those  specified  in  Problem  3.1. 

3.5.2  Suggested  Approach. 

The  electromagnetic  model  can  be  taken  over  directly  from 
Problem  3.4,  assuming  that  a  sufficiently  fine  spatial  grid  has  been 
used.  The  emitting  surfaces  should  run  from  near  the  inlet  (five  cells 
from  the  inlet  boundary)  to  the  corner  of  the  cathode,  and  from  there  to 
the  axis.  The  input  data  for  this  simulation  is  shown  in  Table  3.5. 

3.5.3  Analytical  Solution. 

At  early  time,  the  transient  response  will  be  dominated  by 
plasma  effects  in  the  coaxial  section.  These  results  must  agree  with 
those  from  Problem  3.1;  i.e.,  the  voltage  is  initially  reduced  from  the 
incident  wave  value.  This  effect  is  illustrated  in  the  inlet  voltage 
shown  In  Figure  3.17. 


In  cold-test  (Problem  3.4),  the  diode  presented  an  infinite  im¬ 
pedance.  In  the  present  case,  the  electrons  will  continue  to  flow  across 
the  diode,  since  the  magnetic  field  on  axis  must  vanish  (by  symmetry)  and 
self-insulation  cannot  occur.  Thus,  the  diode  impedance  will  be  reduced. 
An  estimate  of  the  impedance  can  be  obtained  by  assuming  Ch  i  ld-Langmuir 
physics  to  be  applicable.  We  let  the  effective  emission  area  be 


(3.15) 
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Then,  from  Equation  (2.47),  the  total  current  is  given  by 
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Therefore,  the  effective  impedance  must  be  approximately  given  by 
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(d> ) 


where  we  have  used  V  »  d>+.  Using  this  method,  we  estimate  the  impedance 
Z  to  be  30  Q.  This  estimate  can  be  compared  with  the  value  of  29  3 
obtained  from  the  simulation  using  the  expression, 

*3  =  Zi  — -  ,  (3. It  ) 

2d)  -  d> 

where  d>  =  3.2  MV  is  estimated  from  the  simulation.  This  is  good  agree¬ 
ment;  however,  there  is  strong  sensitivity  to  the  choice  of  radius,  R. 


Electron  trajectory  plots  at  two  times  during  the  simulation 
(2.5  and  5.0  nse  )  are  shown  in  Figures  3.18  and  3.19. 
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Table  3.5.  Input  data  for  Problem  3.5. 


title  ‘problem  3.5*  / 

comment  ‘problem  3.5  is  the  same  diode  as  problem  3.4, 
but  with  particles . *  / 
system  2  / 

xlgrid  1  67  2  0.0  25  0.008  0.2  20  0.008  0.1 

20  0 . 002  0.1  / 

x2grid  1  42  2  0.0  20  0.008  0.1  20  0.002  0.1  / 
fields  1  3  800  6.25e-12  0.5  0.5  0.0 

4  1.0  0.29912  0.15022  0.11111  / 
courant  00/ 


conductor  cathode  1  2  22  47  22  47  2  / 
conductor  anode  -1  2  42  67  42  67  2/ 
symmetry  axial  12  2  67  2  / 

particles  emission  null  electrons  0061121 

1.0e+7  6 . 25e-5  2.5e7  0.0  0.0  0.0  0.0  1  7  22  47 

particles  emission  null  electrons  0061121 

1.0e+7  6 . 25e-5  2.5e7  0.0  0.0  0.0  0.0  1  47  22  47 
kinematics  10000/ 
currents  64  1.0  / 
forces  0  0.5  1.0  1.0  / 

voltage  incident  radial  1  twod  01. 012  22  2  42  2/ 
function  incident  0  3  0.0  0.0  1.0e-9  4 . 0e6  1.0  4.0e6  / 
function  radial  5  -11/ 


diagnose  spacing  0  1  courant  01/ 
statistics  200  / 

trajectory  200  1  1  0.0  0.6  -0.15  0.45  / 
observe  1  1  twod  1  0.0  1.0  0.0  1.0  2  2  22  2  42 

twod  1  0.0  1.0  0.0  1.0  2  47  22  47  42 

twod  1  0. 0  1. 0  0. 0  1. 0  1  47  2  67  2  / 

range  200  12  2  22  2  42  1 

1  2  47  22  47  42  1 

1  1  47  2  67  21/ 

display  0  0.0  0.6  -0.15  0.45  / 
output  0  / 
timeout  60  1  / 


start  / 


22  / 
2  / 


**>• 


VERSION:  SEPTEMBER  1393 

T_P7E:  95/ 

SIMULATION:  PROBLEM  3 

.  5 

TIME  HISTORY  PLOT 

E2  COMPONENT 

INTEGRATED  ERCM  (2,221  TO 

(2,421 

Figure  3.17.  Inlet  voltage  vs.  time.  Problem  3 


3.6 


DIODE  WITH  ANODE  FIELD  EMISSION. 


3.6.1 


Problem  Description. 


Problem  3.6  is  a  duplication  of  Problem  3.5,  with  the  addition 
of  ion  field  emission  from  the  radial  portion  of  the  anode.  The  same 
measurements  used  in  the  preceding  simulations  should  be  made,  in  addition 
to  trajectory  plots  of  both  species  (electrons  and  ions). 

3.6.2  Suggested  Approach. 

Including  ions  in  a  diode  simulation  can  cause  problems  due  to 
their  relative  duration.  Since  the  relative  velocities  go  approximately 
as  the  sguare  root  of  the  mass  ratio,  ion  particles  can  remain  in  a 
simulation  much  longer  than  electron  particles.  Thus,  if  ions  and 
electrons  are  treated  eguivalently ,  the  resulting  large  number  of  ion 
particles  can  render  a  simulation  prohibitively  expensive.  This  problem 
has  four  common  remedies: 


(1)  use  ions  only  when  they  impact  physics, 

(2)  increase  the  charge- to-mass  ratio  to  increase  ion 
velocities  artifically, 

(3)  reduce  the  number  of  ions  created  per  cell,  and 

(4)  reduce  the  ion  creation  frequency  (time  step  multiple). 

All  of  the  above  may  be  accomplished  using  options  in  the  PARTICLES 
command.  Input  data  for  this  simulation  is  shown  In  Table  3.6. 

3.6.3  Analytical  Solution. 

At  late  time,  when  many  ions  have  had  time  to  cross  the  gap,  a 
steady  state  with  space-charge- limited  flow  from  both  anode  and  cathode 
can  occur.  The  cathode  will  emit  electrons  and  the  anode  will  emit  ions 
as  if  they  are  both  covered  with  dense  plasma. 
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The  presence  of  ions  in  the  gap  will  allow  higher  electron 
current  densities  and  thus  lower  the  effective  impedance.  The  effective 


current  is 


I  =  tR  (3  ■*  3.  )  , 

e  l 


(3.19) 


where  3^  and  3e  are  the  steady-state  electron  and  ion  currents  in  the 
gap.  The  derivation  which  follows  is  based  on  a  treatment  given  by 
Miller. ^  The  current  densities  are 


3  =  -  en  v 

e  e  e 


3^  =  +en  j  V  .  , 


(3.20) 


where  ne  and  n^  are  the  electron  and  ion  densities  respectively,  e  is 
the  electron  charge,  and  ve  and  v^  are  the  electron  and  ion 
velocities.  We  assume  all  particles  enter  the  gap  with  zero  velocity,  so 
that  the  energies  are  given  by 


Y-V-’J 


1/2  m  v  =  e< 
e  e 


1/2  nruv  f  -  e(  4>0-  1>)  , 


(3.21) 
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where  m^  and  mg  are  the  electron  and  ion  masses,  respectively,  and  J 
is  the  potential.  The  ions  are  (in  this  problem)  assumed  to  be  singly 
charged,  4>0  is  the  potential  at  the  anode,  and  both  anode  and  cathode 
are  assumed  to  behave  as  high-density  plasmas,  so  that 


a  <t>  =  o  , 

X 


(3.22) 


at  both  anode  and  cathode  surfaces.  In  a  manner  similar  to  that  used  for 
pure  electron  flow,  we  write  Poisson's  equation, 


t  R.  6.  Miller,  An  Introduction  to  the  Physics  of  Intense  Charged  Par¬ 
ticle  Beams  Plenum  Publishing,  New  York,  (  1952)  p.  53. 
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where  x  is  the  distance  from  the  cathode.  Assuming  steady-state  flow  for 
both  species,  we  write 

,  3  m  1/2  _  1/  2  _  i/  2 

22<t>  =  —  (  —  )  -  a(  <>  -  0)  1, 

e0  2e 


where 


(3.24) 
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This  equation  can  be  simplified  mathematically  by  substitution  and  a  step 
of  integration  to  obtain 


0  3  !■/  2  l/ 2 

(a  4)  2  =  (  _  )(  _£  )U  +  a(  1-(j<)  ]  -  3  , 

y  9  3 


(3.23) 


where  <1  =  <P/ 4>0,  y  =  x/d  (whete  d  is  the  gap  separation),  and  30  is  the 
electron  current  density  for  such  a  gap  in  the  absence  of  ion  flow. 


The  conditions  of  Equation  (3.22)  at  both  anode  and  cathode 


require 


a  =  1 


3  -  i£  (I®  )  . 
9  3 


(3.26) 


A  second  (numerical)  integration  of  Equation  (3.25)  solves  for  3e/30 
to  yield 
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3  =  1.96  3  . 

e  o 


(3.27) 


and  by  Equations  (3.24-)  and  (3.26), 


3.  =  1.96  3  I  _ 
1  o  '  m  . 
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(3.28) 
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where  30  is  the  current  density  of  electrons  without  ion  flow, 
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Thus,  the  effective  impedance  must  be  approximately  given  by 


1.96  3  tcR  2  (l  +  (  _£  )  ) 

o  v  mi 

,  (3.30) 
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where  Z0  is  the  impedance  for  electron  only  space-charge-limited  flow. 
Therefore,  the  impedance  is  roughly  one-half  the  value  without  ion  flow. 
This  is  a  steady-state  value  which,  for  realistic  ratios  of  m^/me, 

will  not  occur  until  late  in  time.  The  initial  states  of  the  gap  will 

resemble  the  case  with  no  ions;  however,  as  the  gap  approaches  steady- 

state  ion  flow,  the  impedance  will  fall  to  the  final  value. 

Trajectory  plots  for  electrons  and  ions  at  the  end  of  the 

simulation  are  shown  in  Figures  3.20  and  3.21,  respectively.  Since  the 
motion  of  the  ions  is  thus  limited,  there  is  not  much  hope  of  observing 
the  effect  of  Equation  (3.30)  without  runnino  longer  in  time. 
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Table  3.6.  Input  data  for  Problem  3.6. 


title  ‘problem  3.6*  / 

comment  ‘problem  3.6  is  the  same  diode  as  problem  3.5, 
but  with  two  species  of  particles . *  / 
system  2  / 

xlgrid  1  67  2  0.0  25  0.008  0.2  20  0.008  0.1 

10  0.002  0.05  10  0.008  0.05  / 

x2grid  1  42  2  0.0  20  0.008  0.1  10  0.002  0.05  10  0.008  0.05  / 

fields  1  3  800  6.25e-12  0.5  0.5  0.0 

4  1.0  0.29912  0.15022  0.11111  / 
courant  00/ 

conductor  cathode  1  2  22  47  22  47  2  / 
conductor  anode  -1  2  42  67  42  67  2  / 
symmetry  axial  12  2  67  2  / 

particles  emission  null  electrons  0061121 

1.0e+7  6 . 25e-5  2.5e7  0.0  0.0  0.0  0.0  1  7  22  47  22  / 

particles  emission  null  electrons  0061121 

1.0e+7  6 . 25e-5  2.5e7  0.0  0.0  0.0  0.0  1  47  22  47  2  / 

particles  emission  null  ions  1111121 

1.0e+7  6 . 25e-5  2 . 5e7  0.0  0.0  0.0  0.0  -1  67  42  67  2  / 

kinematics  10000/ 
currents  64  1.0  / 
forces  0  0.5  1.0  1.0  / 

voltage  incident  radial  1  twod  01. 012  22  2  42  2/ 
function  incident  0  3  0.0  0.0  1.0e-9  4.0e6  1.0  4.0e6  / 
function  radial  5  -11/ 

diagnose  spacing  0  1  courant  01/ 
statistics  200  / 

trajectory  200  1  1  0.0  0.6  -0.15  0.45  / 
observe  1  1  twod  1  0.0  1.0  0.0  1.0  2  2  22  2  42 

twod  1  0.0  1.0  0.0  1.0  2  47  22  47  42 

twod  10. 01. 00. 01. 01  47  2  67  2/ 

range  200  12  2  22  2  42  1 

1  2  47  22  47  42  1 

1  1  47  2  67  21/ 

display  0  0.0  0.6  -0.15  0.45  / 
output  0  / 
timeout  60  1  / 
start  / 
stop  / 
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3.7 


COLD-TEST  OF  SHORTED  DIODE. 


3.7.1  Problem  Description. 

In  this  simulation,  a  short  is  to  be  placed  between  anode  and 
cathode  two  cells  above  the  axis  of  symmetry,  and  the  field  emission  is  to 
be  deleted.  Measurements  consisting  of  range  plots  and  time  histories 
taken  at  the  inlet,  the  mismatch,  and  near  the  short,  should  be  made. 

3.7.2  Suggested  Approach. 

Maintaining  the  previous  electromagnetic  model,  simply  delete 
the  axial  symmetry  boundary  condition  and  insert  a  conductor.  The  input 
data  is  shown  in  Table  3.7. 

3.7.3  Analytical  Solution. 

The  results  derived  for  Problem  3.2  are  clearly  applicable 
here.  The  multiple  reflections  described  in  Equation  (9)  are  particularly 
pronounced  due  to  the  disparity  in  impedances,  i.e.,  li  *  “  and  Z3  =  0. 
Voltage  plots  at  inlet,  midpoint,  and  near  axis  are  shown  in  Figures  3.22, 
3.23,  and  3.24,  respectively. 
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Table  3.7.  Input  data  for  Problem  3.7. 


title  *problem  3.7*  / 

comment  *problem  3.7  is  a  cold  test  of  problem  3.6 
with  a  short*  / 
system  2  / 

xlgrid  1  67  2  0.0  25  0.008  0.2  20  0.008  0.1 

20  0.002  0.1  / 

x2grid  1  42  2  0.0  20  0.008  0.1  20  0.002  0.1  / 
fields  1  3  800  6.25e-12  0.5  0.5  0.0 

4  1.0  0.29912  0.15022  0.11111  / 
courant  00/ 

conductor  cathode  1  2  22  47  22  47  4  / 

conductor  short  1  47  4  67  4  / 

conductor  anode  -1  2  42  67  42  67  4/ 

symmetry  axial  12  2  67  2  / 

z  particles  emission  null  electrons  0061121 

1.0e+7  6 . 25e-5  2.5e7  0.0  0.0  0.0  0.0  1  7  22  47  22  / 

z  particles  emission  null  electrons  0061121 

1.0e+7  6 . 25e-5  2 . 5e7  0.0  0.0  0.0  0.0  1  47  22  47  4  / 

kinematics  10000/ 
currents  64  1.0  / 
forces  0  0.5  1.0  1.0  / 

voltage  incident  radial  1  twod  01. 012  22  2  42  2/ 
function  incident  0  3  0.0  0.0  l.Qe-9  4.0e6  1.0  4.0e6  / 
function  radial  5  -11/ 

diagnose  spacing  0  1  courant  01/ 
statistics  200  / 

trajectory  200  1  1  0.0  0.6  -0.15  0.45  / 
observe  1  1  twod  1  0.0  1.0  0.0  1.0  2  2  22  2  42 

twod  1  0.0  1.0  0.0  1.0  2  47  22  47  42 

twod  10. 01. 00. 01. 01  47  5  67  5/ 

range  200  12  2  22  2  42  1 

1  2  47  22  47  42  1 

1  1  47  5  67  51/ 

display  0  0.0  0.6  -0.15  0.45  / 
output  0  / 
timeout  60  1  / 
start  / 
stop  / 
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3.8 


SHORTED  DIODE  WITH  FIELD  EMISSION. 


3.8.1  Problem  Description. 

Problem  3.S  uses  the  shorted  diode  of  Problem  3.7,  but  with 

cathode  (electron)  field  emission  reinstated  on  both  conformal  surfaces. 

3.3.2  Suggested  Approach. 

A  potential  difficulty  can  result  from  electron  trajectory 

calculations  due  to  the  large  magnetic  fields  near  axis.  The  standard 
kinematics  algorithm  will  allow  time  steps  consistent  with  up  to  it/4 

radians  in  orbit.  Beyond  this,  the  algorithm  becomes  inaccurate,  and 

artificial  particle  motion  will  result.^-  The  onset  of  this  difficulty 
can  be  estimated  using  previous  results  from  Problems  2.2  and  3.7. 

Input  data  for  this  simulation  is  shown  in  Table  3.8. 

3.8.3  Analytical  Solution. 

Current  through  the  short  produces  a  magnetic  field  which  varies 
inversely  with  radius.  Thus,  electrons  tend  to  move  in  circular  orbits, 
and  tend  to  move  toward  the  short  due  to  ExB  drift  (see  Problem  2.3). 
However,  the  short  impedance  is  dominant,  so  that  field  emission  affects 
the  transient  but  not  the  steady  state. 

Figure  3.25  illustrates  the  voltage  at  inlet.  Electron 
trajectory  plots  at  2.5  and  5.0  nsec  are  shown  in  Figures  3.26  and  3.27, 
respectively,  where  the  tendency  of  electrons  to  orbit  is  clearly 
demonstrated . 


t  B.  Gopien.  3.  Brandenburg,  and  R.  Worl,  "Particle  Subcvcling  in  Pulsed 
Power  Simulations."  MRC/WDC-R- 125.  April  1987. 
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Table  3.S.  Input  data  for  Problem  3.S. 


title  ^problem  3.8*  / 

comment  *problem  3.8  is  a  simulation  of 

a  diode  like  the  one  in  problem  3.5, 
but  with  a  short  near  the  axis  of  symmetry*  / 
system  2  / 

xlgrid  1  67  2  0.0  25  0.008  0.2  20  0.008  0.1 

20  0.002  0.1  / 

x2grid  1  42  2  0.0  20  0.008  0.1  20  0.002  0.1  / 
fields  1  3  800  6.25e-12  0.5  0.5  0.0 

4  1.0  0.29912  0.15022  0.11111  / 
courant  00/ 

conductor  cathode  1  2  22  47  22  47  4  / 

conductor  short  1  47  4  67  4  / 

conductor  anode  -1  2  42  67  42  67  4/ 

symmetry  axial  12  2  67  2  / 

particles  emission  null  electrons  0061121 

1.0e+7  6 . 25e-5  2.5e7  0.0  0.0  0.0  0.0  1  7  22  47  22  / 

particles  emission  null  electrons  0061121 

1.0e+7  6 . 25e-5  2.5e7  0.0  0.0  0.0  0.0  1  47  22  47  4 

kinematics  10000/ 
currents  64  1.0  / 
forces  0  0.5  1.0  1.0  / 

voltage  incident  radial  1  twod  01. 012  22  2  42  2/ 
function  incident  0  3  0.0  0.0  1.0e-9  4.0e6  1.0  4.0e6  / 
function  radial  5  -11/ 

diagnose  spacing  0  1  courant  01/ 
statistics  200  / 

trajectory  200  1  1  0.0  0.6  -0.15  0.45  / 
observe  1  1  twod  1  0.0  1.0  0.0  1.0  2  2  22  2  42 

twod  1  0.0  1.0  0.0  1.0  2  47  22  47  42 

twod  10. 01. 00. 01. 01  47  5  67  5/ 

range  200  12  2  22  2  42  1 

1  2  47  22  47  42  1 

1  1  47  5  67  51/ 

display  0  0.0  0.6  -0.15  0.45  / 
output  0  / 
timeout  60  1  / 
start  / 
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iqure  3.26.  Trajectory  plot  at  2.5  nsec.  Problem  3.8. 


rRG I C  VERSION:  SEPTEMBER  '.933  ERIE:  SS/G3/:3 
simulation:  problem  3.3 


igure  3.27.  Trajectory  plot  at  5.0  nsec.  Problem  3.8. 
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